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Abstract 


A  series  of  tables  and  charts  is  presented  from  which  the  atmospheric 
transmittance  between  any  two  points  in  the  terrestrial  atmosphere  can  be  deter¬ 
mined.  This  material  is  based  on  a  set  of  five  atmospheric  models  ranging  from 
tropical  to  arctic  and  two  aerosol  models.  A  selected  set  of  laser  frequencies 
has  been  defined  for  which  monochromatic  transmittance  values  have  been  given. 
For  low  resolution  transmittance  prediction,  a  series  of  charts  has  been  drawn 
providing  the  capability  for  predicting  transmittance  at  a  resolution  of  20  wave- 
numbers.  Separate  sections  are  included  on  scattered  solar  radiation,  infrared 
emission,  refractive  effects,  and  attenuation  by  cloud  and  fog. 
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Optical  Properties  of  the  Atmosphere 

(Revised) 

1.  INTRODUCTION 

This  report  is  intended  to  provide  data  from  which  the  attenuation  or  emission 
of  radiation  along  a  path  in  a  given  model  atmosphere  can  be  determined.  The 
spectral  region  covered  is  from  0.25  to  25  micrometers.  We  have  not  concerned 
ourselves  with  the  atmosphere  above  100  km,  thereby  neglecting  less  than  one- 
millionth  of  the  total  atmospheric  mass. 

We  have  defined  five  model  atmospheres  each  for  temperature,  pressure  and 
absorbing  gas  concentrations.  In  addition  to  these,  two  aerosol  models  have  been 
defined.  Thus,  one  effectively  has  10  models  from  which  to  choose  In  order  to  use 
the  data  as  presented  without  interpolation.  However,  Interpolation  schemes  are 
described  in  order  to  allow  the  user  to  obtain  optical  data  for  intermediate  condi¬ 
tions. 

As  It  is  not  possible  to  cover  all  conceivable  atmospheric  models,  neither  Is 
It  possible  to  provide  data  in  a  short  report  for  any  spectral  resolution.  Thus,  we 
are  forced  to  limit  the  high  resolution  calculations  to  a  few  specific  wavelengths. 

A  list  of  laser  emission  wavelengths  was  chosen  for  this  purpose.  We  attempted  to 
pick  those  lasers  most  commonly  used  In  practical  applications  and  as  research 
tools.  In  addition  to  the  laser  beam  transmittance  table,  we  have  provided  in 
graphical  form  a  calculation  scheme  for  determining  atmospheric  absorption  and 
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mattering  at  about  a  110  wavenumber  resolution  across  the  entire  apeotral 

ration. 

The  models  defined  in  Meotion  3  together  with  the  low  resolution  absorption 
calculation  technique*  presented  in  Section  8.  were  used  to  compute  the  radiation 
emitted  to  apace  and  the  radiation  aa  observed  looking  up  at  the  sky  from  the  sur¬ 
face  (aky  background)  for  wavelengths  longer  than  8pm  (Section  6).  The  radiation 
to  apace  and  aky  backgrounds  at  shorter  wavelengths  ts  a  complicated  function  of 
the  solar  aenith  angle  and  the  angle  of  observation,  so  only  a  few  specific  examples 
have  been  given. 

The  effect  of  the  refraotive  index  on  the  geometrical  path  length  and  the  effect 
of  cloud  and  fog  on  the  attenuation  of  radiation  are  discussed  in  Sections  7  and  8, 
respectively. 

8,  MimKI.  M  MOSPIIKHtlS  USKI)  IN  OPTICAL  CALCULATIONS 
3.)  Slundard  MuiM  tlmoNphere* 

A  aeries  of  five  model  atmospheres  is  presented  in  Table  1.  The  pressure, 
temperature,  density,  water  vapor  and  ozone  concentrations  are  those  of  the  U.S. 
Standard  Atmosphere  of  1863  and  the  Supplemental  Atmospheres  as  compiled  in 
the  Handbook  of  Geophysics  and  Space  Environment  (Valley,  1965).  The  water 
vapor  densities  above  11  km  have  been  taken  from  Sissenwlne  et  al  (1968).  A  dis¬ 
cussion  of  units  and  conversion  can  be  found  in  the  Appendix. 


3.3  Concentration*  of  Uniformly  Mixed  (iu«o* 

Table  3  1b  intended  to  provide  the  reader  with  the  data  for  determining  the 
concentrations  of  all  uniformly  mixed  atmospheric  gases  in  any  horizontal  or  verti¬ 
cal  path.  Column  3  gives  the  concentrations  in  parts  per  million  by  volume.  The 
fourth  column  gives  the  totsl  amount  of  gas  in  a  vertical  path  through  the  entire 
atmosphere.  The  fifth  column  lists  the  number  of  (cm-atm)STp  of  gaB  in  a  one 
kilometer  path  at  sea  level.  The  number  of  (cm-atm)g^.p  in  a  horizontal  path  at 
any  other  altitude  is  given  by  Eq.  (1). 


(cm-atm)g.pp 


(1) 


where  L  is  the  path  length  in  km,  and  PQ  and  TQ  refer  to  STP  (PQ  *  1  atm  and 
T0  »  273K). 


Table  1.  Model  Atmosphere*  Used  as  a  Basts  of  the  Computation 
of  Atmospheric  Optical  Properties 


TROPICAL 

1 

Ht. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 

(g/m* 

Water  Vapor 
(g/m3) 

Oaone 

(g/m3) 

0 

1.  013E+03 

300.  0 

1.  167E+03 

T.  9E+01 

5.  6E-05 

1 

9.  040E+02 

294.  0 

1.  064E+03 

1.  3E+01 

5.  6E-05 

2 

8.  050E+02 

288.  0 

9.  689E+02 

9.  3E+00 

5.4E-05 

3 

7.  150E+02 

284.  0 

8.  756E+02 

4.7E+00 

5.  IE-05 

4 

6.  330E+02 

277.  0 

7.  951E+02 

2.2E+00 

4.  7E-05 

5 

5.  590E+02 

270.  0 

7.  199E+02 

1.  5E+00 

4.  5E-05 

6 

4.  920E+02 

264.  0 

6.  501E+02 

8.  5E-01 

4.  3E-05 

7 

4.  320E+02 

257.  0 

5.  855E+02 

4.  7E-01 

4.  IE -05 

8 

3.780E+02 

250.  0 

5.  258E+02 

2.  5E-01 

3.  9E-05 

9 

3. 290E+02 

244.  0 

4.  708E+02 

1.  2E-01 

3.  9E-05 

10 

2. 860E+02 

237.  0 

4.  2  02E+02 

5.  0E-02 

3.  9E-05 

11 

2. 470E+02 

230.  0 

3.  740E+02 

1.  7E-02 

4.  IE-05 

12 

2.  130E+02 

224.  0 

3.  316E+02 

6.  0E-03 

4.  3E-05 

13 

1.  820E+02 

217.  0 

2.  929E+02 

1.  8E-03 

4.  5E-05 

14 

1.  560E+02 

210.  0 

2.  578E+02 

1.  0E-03 

4.  5E-05 

15 

1.  320E+02 

204.  0 

2.  260E+02 

7.  6E-04 

4.  7E-05 

16 

1.  110E+02 

197. 0 

1.  972E+02 

6.  4E-04 

4.  7E-05 

17 

9.  370E+01 

195.  0 

1.  676E+02 

5.  6E-04 

6.  9E-05 

16 

7.  890E+01 

199.  0 

1.  382E+02 

5.  0E-04 

9.  0E-05 

19 

6,  660E+01 

203.  0 

1.  145E+02 

4.  9E-04 

1.  4E-04 

20 

5.  650E+01 

207.  0 

9.  515E+01 

4.  5E-04 

1.  9E-04 

21 

4.  800E+01 

211.  0 

7. 938E+01 

5.  IE-04 

2.  4E-04 

22 

4.  090E+01 

215.  0 

6.  645E+01 

5.  IE -04 

2.  8E-04 

23 

3.  500E+01 

217.  0 

5.  618E+01 

5.  4E-04 

3.  2E-04 

24 

3.  OOOE+Ol 

219.  0 

4.  763E+01 

6.  OE-04 

3.  4E-04 

25 

2.  570E+01 

221.  0 

4. 045E+O1 

6.  7E-04 

3.  4E-04 

30 

1.  220E+01 

232.  C 

1.  831E+01 

3.  6E-04 

2.  4E-04 

35 

6.  OOOE+OO 

243.  0 

8.  600E+00 

1.  IE -04 

9.  2E-05 

40 

3.  050E+00 

254.  0 

4.  181E+00 

4.  3E-05 

4.  IE-05 

45 

1.  590E+00 

265.  0 

2.  097E+00 

1.  9E-05 

1.  3E-05 

50 

8.  540E-01 

270.  0 

1.  101E+00 

6.  3E-06 

4.  3E-06 

70 

5.  790E-02 

219.0 

9.  2 10E  -02 

1.4E-07 

8.  6E-08 

100 

3.  000E-04 

210.  0 

5.  000E-04 

1.  0E-09 

4,  3E-11 

4 
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Preaaure 

(mb) 


1.013E+03 
9.  020E+02 
8.  020E+02 

7.  100E+02 

6.  280E+02 

5.  54QE+02 

4.  870E+02 
4. 260E+02 
3. 720E+02 

3.  240E+02 

2.  810E+02 
2. 430E+02 
2. 090E+02 

1.  7  90E+02 

1.  530E+02 
1.  300E+02 
l.  110E+02 
9.  500E+01 

8.  120E+01 

6.  950E+01 

5.  950E+01 

5.  100E+01 

4.  370E+01 

3.  760E+01 

3.  220E+01 

2.  770E+01 
1.  320E+01 

6.  520E+00 

3.  330E+00 
1.  760E+00 

9.  510E-01 
6.  7 10E-02 

3.  000E-04 


MIDLATITUDE  SUMMER 


Tamp. 

(°K) 

294.  C 
290,  0 
288.  0 
279.  0 
273.  0 
267.  0 
261.  0 
255.0 
248.  0 
242.  0 
235.  0 
229.  0 
222.  0 
216.  0 
216.  0 
216.  0 
216.  0 

215.  0 

216.  0 
217. 0 
218.  0 

219.  0 

220.  0 
222.  0 

223.  0 

224.  0 
234.  0 
243.  0 
258.  0 
270.  0 
276.  0 
218.  0 
210.  0 


Dfnalty 

_ ta/mfry 

1.  181E+03 
1.  080E+03 
9.  757E+02 
8.  846E+02 

7.  998E+02 

7.  211E+02 
6.  487E+02 

5.  830E+02 
5.  225E+02 

4.  669E+02 

4.  159E+02 
3.  693E+02 
3.  269E+02 

2.  882E+02 
2.  464E+02 
2.  104E+02 
1.  797E+02 
1.  535E+02 
1.  305E+02 
1.  110E+02 
9.  453E+01 

8.  056E+01 

6.  872E+01 

5.  867E+01 

5.  014E+01 
4. 288E+01 
1.  322E+01 

6.  5  19Ef00 
3. 330E+00 
1.  757E+00 

9.  512E-01 
0.  706E-02 
5.  000E-04 


XrJrr 

1.4E+01 
9.  3E+00 

5.  9E+00 
3.  3E+00 
1.  9E+00 

1.  0E+00 

6.  IE-01 
3.7E-01 

2.  IE-01 
1.2E-01 
6.  4E-02 
2.  2E-02 

6.  0E-03 
1.  8E-03 
1.  0E-03 

7.  6E-04 
6.  4E-04 
5.  6E-04 
5.  0E-04 

4.  9E-04 

4.  5E-04 

5.  IE-04 

5.  IE -04 

5,  4E-04 

6.  OF -04 

6.  7E-04 

3.  6E-04 
1.  IE -04 

4.  3E-05 
1.  9E-05 
6.  3E-06 
1.4E-07 
1.  0E-09 


Oaong 

<g/m3) 

6.  0E-05 
C.0E-05 
6.  QE-05 
6.2E-05 
6.  4E-05 
6.  6E-05 

6.  9E-05 

7.  SE-05 

7.  9E-05 

8.  6E-05 

9.  0E-05 
1.  IE -04 
1.2E-04 
1.  5E-04 
1.  8E-04 

1.  9E-04 

2.  IE -04 
2.4E-04 

2.  8E-04 

3.  2E-04 
3.  4E-04 
3.  6E-04 
3.  6E-04 
3.  4E-04 
3.2E-04 

3.  OE-04 
2.  0E-04 
9.  2E-05 

4.  IE-05 
1.  3E-05 
4,  3E-06 
8.  6E-08 
4.  3E-11 
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Table  1  (Contd).  Model  Atmosphere!  Used  ie  a  Beale  of  the  Computation 
of  Atmoapherlc  Optical  Properties 


MIDLATITUDE  WINTER 

Ht 

(km) 

Pressure 

(mb) 

Temp. 

(OK) 

Density 
(g/m  3) 

Water  Vapor 
(g/m3) 

Ozone 
(g/m  3) 

0 

1.  018E+03 

272.2 

1.  301E+03 

3.  5E+00 

6.  0E-05 

1 

8.  973E+02 

268.  7 

1.  162E+03 

2.  5E+00 

5.4E-05 

2 

7.  897E+02 

265.2 

1.  037E+03 

1.  8E+00 

4.  9E-05 

3 

6.  938E+02 

261.  7 

9.  230E+02 

1.2E+00 

4.  9E-05 

4 

6.  081E+02 

255.7 

8.282E+02 

6.  6E •  0 1 

4.  9E-05 

5 

5.  313E+02 

249.  7 

7,41 1E+02 

3.  8E-01 

5.  8E-05 

6 

4.  627E+02 

243.  7 

6.  614E+02 

2.  IE-01 

6.4E-05 

7 

4.  016E+02 

237.7 

5.  886E+02 

8.  5E-02 

7.7E-05 

8 

3.  473E+02 

231.  7 

S.  222E+02 

3.  5E-02 

9.  0E-05 

9 

2.  992E+02 

225.  7 

4.  619E+02 

1.  6E-02 

1.2E-04 

10 

2.  568E+02 

219.7 

4.  072E+02 

7.  5E-03 

1.  6E-04 

11 

2.  199E+02 

219.  2 

3.496E+02 

6.  9E-03 

2.  IE -04 

12 

1.  882E+C2 

218.7 

2. 999E+02 

6.  0E-03 

2.  6E-04 

13 

1.  610E+02 

2  18.2 

2. 572E+02 

1.  8E-03 

3.  0E-04 

14 

1.  378E+02 

217.  7 

2. 206E+02 

1.  0E-03 

3.  2E-04 

15 

1.  178E+02 

217.2 

1.  890E+02 

7.  6E-04 

3.4E-04 

16 

1.  007E+02 

216.  7 

1.  620E+02 

6.4E-04 

3.  6E-04 

17 

8.  610E+01 

216.  2 

1.  388E+02 

5.  6E-04 

3.  9E-04 

18 

7. 350E+01 

215.  7 

1.  188E+02 

5.  0E-04 

4.  IE-04 

19 

6.2  80E+01 

215.  2 

1.  0 17E+02 

4.  9E-04 

4.  3E-04 

20 

5.  370E+01 

215.2 

8.  690E+01 

4.  5E-04 

4.  5E-04 

21 

4.  580E+01 

215.  2 

7.421E+01 

5.  IE -04 

4.  3E-04 

22 

3.  910E+01 

215.  2 

6.  338E+01 

5.  IE-04 

4.  3E-04 

23 

3.  340E+01 

215.2 

5.415E+01 

5.  4E-04 

3.  9E-04 

24 

2. 860E+01 

2  15.  2 

4.  624E+01 

6.  0E-04 

3.  6E-04 

25 

2 . 430E+0 1 

215.  2 

3.  950E+01 

6.  7E-04 

3.  4E-04 

30 

1.  110E+01 

217.  4 

1.  783E+01 

3.  6E-04 

1.  9E-04 

35 

5. 180E+00 

227.  8 

7.  924E+00 

1.  IE -04 

9.2E-05 

40 

2. 530E+00 

243.  2 

3,  625E+00 

4,  3E-05 

4,  IE-05 

45 

1.  290E+00 

258.  5 

1.  741E+00 

1.  9E-05 

1.  3E-05 

50 

6.  820E-01 

265.7 

8.  954E-01 

6.  3E-06 

4.  3E-06 

70 

4.  670E-02 

230.7 

7.  051E-02 

1.  4E-07 

8.  6E-08 

100 

3.  000E-04 

210.2 

5.  000E-04 

1.  0E-09 

4.3E-11 
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Table  1  (Contd).  Model  Atmoepherea  Used  as  a  Baala  of  the  Computation 
of  Atmoepherio  Opttoal  Propertiea 


SUBARCTIC  SUMMER 

Ht. 

(tan) 

Preaeure 

(mb) 

Temp. 

(°K) 

BSElSflHI 

Water  Vapor 
<g/m*> 

Ozone 

(g/m3) 

0 

1.  010E+03 

287.0 

1.  220E+03 

8.  1E+00 

4.  9E-05 

1 

8.  960E+02 

282.0 

1.  110E+03 

6.  0E+00 

5.4E-05 

2 

7.  929E+02 

276.0 

9.  97 1E+02 

4.  2E+00 

5.6E-05 

3 

7.  OOOE+02 

271.0 

8.  985E+02 

2.7E+00 

5.  8E-05 

4 

6.  160E+02 

266.0 

8.  077E+02 

1.  7E+00 

6.  0E-05 

5 

5.  4 10E+02 

260.  0 

7. 244E+02 

1.  0E+00 

6.4E-05 

6 

4.  730E+02 

253.  0 

6.  519E+02 

5.4E-01 

7.  IE-05 

7 

4.  130E+02 

246.  0 

5.  849E+02 

2.  9E-01 

7.  5E-05 

8 

3.  590E+02 

239.0 

5.231E+02 

1.  3E-02 

7.  9E-05 

9 

3.  107E+02 

232.0 

4.  663E+02 

4.2E-02 

1.  IE -04 

10 

2.  677E+02 

225.0 

4.  142E+02 

1.  5E-02 

1.  3E-04 

11 

2.  300E+02 

225.0 

3.  559E+02 

9.  4E-03 

1.  8E-04 

12 

1.  977E+02 

225.0 

3.  059E+02 

6.  0E  -03 

2.  IE -04 

13 

1.  700E+02 

225.  0 

2.  630E+02 

1.  8E-03 

2.  6E-04 

14 

1.  460E+02 

225.  0 

2.  260E+02 

1.  0E  -03 

2.  8E-04 

15 

1.  250E+02 

225.  0 

1.  943E+02 

7.  6E-04 

3.  2E-04 

16 

1.  080E+02 

225.0 

1.67  1E+02 

6.  4E-04 

3.4E-04 

17 

9.  280E+01 

225.  0 

1.  436E+02 

5.  6E-04 

3.  9E-04 

18 

7.  980E+01 

225.  0 

1.  235E+02 

5.  0E-04 

4.  IE -04 

19 

6.  860E+01 

225.  0 

1.  062E+02 

4.  9E-04 

4.  IE -04 

20 

5.  800E+01 

225.  0 

9.  128E+01 

4.  5E-04 

3.  9E-04 

21 

5.  070E+01 

225.  0 

7. 849E+01 

5.  IE -04 

3.  6E-04 

22 

4.  360E+01 

225.  0 

6.  750E+01 

5.  IE -04 

3.2E-04 

23 

3. 750E+01 

225.  0 

5.  805E+01 

5.4E-04 

3.  0E-04 

24 

3.  227E+01 

226.  0 

4. 963E+01 

6.  0E-04 

2.  8E-04 

25 

2.  780E+01 

228.  0 

4. 247E+01 

6.  '"E-04 

2.  6E-04 

30 

1.  340E+01 

235.0 

1.  338E+01 

3.  6E-04 

1.4E-04 

35 

6.  610E+00 

247.0 

6.  614E+00 

1.  IE -04 

9.  2E-05 

40 

3.  400E+00 

262.  0 

3.  404E+00 

4.  3E-05 

4.  IE-05 

45 

1.  810E+00 

274.0 

1.  817E+00 

1.  9E-05 

1.  3E-05 

50 

9.  870E-01 

277.0 

9.  868E-01 

6.  3E-06 

4.  3E-06 

70 

7.070E-02 

216.0 

7.07  IE-02 

1.4E-07 

8.  6E-08 

100 

3.  OOOE-04 

210.0 

5.  000E-04 

1.  0E-09 

4.  3E-11 

Table  1  (Contd) .  Model  Atmospheres  Used  as  a  Basis  of  the  Computation 
of  Atmospheric  Optical  Properties 


SUBARCTIC  WINTER 

Ht. 

Pressure 

Temp. 

Density 

Water  Vapor 

Ozone 

(km) 

(mb) 

(°K) 

(g/m3) 

(g/m3) 

(g/m3) 

0 

1.  013E+03 

257.  1 

1.  372E+03 

1.  2E+00 

4.  IE-05 

1 

8.  878E+02 

259.  1 

1.  193E+03 

1.  2E+00 

4.  IE-05 

2 

7.  775E+02 

255.  9 

1.  058E+03 

P.4E-01 

4.  IE-05 

3 

6.  798E+02 

252.7 

9.  366E+02 

6.  8E-01 

4.  3E-05 

4 

5.  932E+02 

247.7 

8.  339E+02 

4. IE-01 

4.  5E-05 

5 

5.  158E+02 

240.  9 

7.  457E+02 

2. 0E-01 

4.  7E-05 

6 

4.  467E+02 

234.  1 

6.  646E+02 

9.  8E-02 

4.  9E-05 

7 

3.  853E+02 

227.  3 

5.  904E+02 

5.4E-02 

7.  IE-05 

8 

3.  308E+02 

220.  6 

5.226E+02 

1.  IE -02 

9.  0E-05 

8 

2. 829E+02 

217.2 

4.  538E+02 

8.4E-03 

1.  6F-04 

10 

2. 418E+02 

217.2 

3.  879E+02 

5.  5E-03 

2.  4E-04 

11 

2.  067E+02 

217.2 

3.  315E+02 

3.  8E-03 

3.2E-04 

12 

1.  766F+02 

217.2 

2.  834E+02 

2.  6E-03 

4.  3E-04 

13 

1.  510E+02 

217.  2 

2.  422E+02 

1.  8E-03 

4.  7E-04 

14 

1.  291E+02 

217.2 

2.  071E+02 

1.  0E-03 

4.  9E-04 

15 

1.  103E+02 

217.2 

1.  770E+02 

7.  6E-04 

5.  6E-04 

16 

9.431E+01 

216.6 

1.  517E+02 

6.4E-04 

6.2E-04 

17 

8. 058E+01 

216.0 

1.  300E+02 

5.  6E-04 

6.2E-04 

18 

6.  882E+01 

215.4 

1.  113E+02 

5.  0E-04 

6.2E-04 

19 

5.  875E+01 

214.8 

9.  529E+01 

4.  9E-04 

6.  0E-04 

20 

5.  014E+01 

214.  1 

8.  155E+01 

4.  5E-04 

5.  6E-04 

21 

4.  277E+01 

213.  6 

6,  976E+01 

5.  IE -04 

5.  IE -04 

22 

3.  647E+01 

213.0 

5.  966E+01 

5.  IE -04 

4.7E-04 

23 

3.  109E+01 

212.4 

5.  100E+01 

5.4E-04 

4.  3E-04 

24 

2.  649E+01 

211.  8 

4.  358E+01 

6.  0E-04 

3.  6E-04 

25 

2.  256E+01 

211.2 

3.  722E+01 

6.  7E-04 

3.2E-04 

30 

1.020E+01 

216.  0 

1.  645E+01 

3.  6E-04 

1.  5E-04 

35 

4. 701E+00 

222.2 

7.  368E+00 

1.  IE -04 

9.2E-05 

40 

2.  243E+00 

234.7 

3.  330E+00 

4.3E-05 

4.  IE -05 

45 

1.  113E+00 

247.0 

1.  569E+00 

1.  9E-05 

1.  3E-05 

50 

5.719E-01 

259.3 

7.  682E-01 

6.  3E-06 

i.  3E-06 

70 

4.  016E-02 

245.7 

5.  695E-02 

1.4E-07 

8.  6E-08 

100 

3.  000E-04 

210.0 

5.  Q00E-04 

1.  0E-09 

4.3E-11 

Table  2.  Concentrations  of  Uniformly  Mixed  Gases 


Constitu* 

ent 

Mole¬ 

cular 

Wt. 

ppm  by 
Vol. 

(cm-atm)gTp 

in  vertical  path 
from  sea  level 

(cm-atm)STpykm 

in  horizontal 
path  at  sea  level 

gm  cm"2/mb 

Ref. 

Air 

28.97 

106 

8xl05 

10® 

1.  02 

* 

co2 

44 

330 

264 

33 

5.  JlxlO-4 

t 

n2o 

44 

0.  28 

0.22 

0.028 

4.  34xl0'7 

+ 

CO 

28 

0.  075 

0.  06 

0. 0075 

7.  39xl0-8 

0 

ch4 

16 

1.  6 

1.28 

0.  16 

9.  01xl0‘7 

** 

°2 

32 

2.  095xl05 

1.  68x10® 

2. 095xl04 

0.236 

* 

The  sixth  column  provides  values  from  which  the  amount  of  gas  between  any 
two  pressure  levels  in  the  atmosphere  can  be  determined  by  multiplication  of  the 
tabulated  values  by  the  pressure  increment  in  millibars  (Ap)  times  the  secant  of 
the  zenith  angle  (see  Appendix). 

2.3  Aerosol  Distributions 

The  two  aerosol  models  (Table  3)  describe  a  "clear"  and  a  "hazy"  atmosphere 
corresponding  to  a  visibility  of  23  and  5  km  respectively  at  ground  level.  The 
aerosol  size  distribution  function  for  both  models  is  the  Bame  at  all  altitudes  and 
is  similar  to  the  one  suggested  by  Deirmendjian  (1963)  for  continental  haze 
(Figure  1),  whereby  the  relationship  between  the  number  density  (n(r)j  and  the 
particle  radius  is  given  by: 

-  C.-r'4  for  0.  1  pm  <  r  <  10  pm 
Ar  i 

■  C  j*  104  for  0.  02  pm  <  r  <  0.  1  pm  (2) 

■  o  for  r  <  0. 02  and  r  >  10. 0. 

It  differs  from  Delrmendjian's  model  "C"  in  that  the  Urge  particle  cutoff  has  been 
extended  from  S  pm  to  10  pm.  This  distribution  function  is  normalised  with 
Ar  ■  1  pm  and  ■  0.  883  x  10*3  ,  (Figure  1),  that  is.  the  Integral  over  the  whole 
else  range  becomes  1.  One  need  only  multiply  the  normalised  distribution  function 

•  Valley  (1985),  *  Fink  et  al  (1964>*Btrkland  and  Shaw  (1959),  °Shaw  (1968), 
“Shaw  (1989). 


Table  3.  Aerosol  Models  -  Vertical  Distributions  for  a  "Clear" 
and  "Hazy"  Atmosphere 


PAR  TIC  LE  DENSI T  Y  , 

N  (PARTICLES  PER  cmd) 

Altitude 

23  km  Visibility 

5  km  Visibility 

(km) 

Clear 

Hazy 

0 

2.  828E+03 

1.  378E+04 

1 

1.  244E+03 

5.  030E+03 

2 

5.  37  1E+02 

1.844E+03 

3 

2.  256E+02 

6.  731E+02 

4 

1.  192E+02 

2.  453E+02 

5 

8.  987E+01 

8.  987E+01 

6 

6.  337E+01 

6.  337E+01 

7 

5.  890E+01 

5.  890E+01 

8 

6.  069E+01 

6.  069E+01 

9 

5.  818E+01 

5.  818E+01 

10 

5.  675E+01 

5.  675E+01 

11 

5.  317E+01 

5.  317E+01 

12 

5.  585E+01 

5.  585E+01 

13 

5.  156E+01 

5.  156E+01 

14 

5.  048E+01 

5.  048E+01 

15 

4.  744E+01 

4.  744E+01 

16 

4.  511E+01 

4.  511E+01 

17 

4.  458E+01 

4.  453E+01 

18 

4.  314E+01 

4.  314E+01 

19 

3.  6342+01 

3.  634E+01 

20 

2.  667E+01 

2.  667E+01 

21 

1.  933E+01 

1.  933E+01 

22 

1.  455E+01 

1.  455E+01 

23 

1.  118E+01 

1.  113E+01 

24 

8.  826E+00 

8.  826E+00 

25 

7.  429E+00 

7.  429E+00 

30 

2.  238E+00 

2.  238E+00 

35 

5.  890E-01 

5.  890E-01 

40 

1.  550E-01 

1.  550E-01 

45 

A.  082E-02 

4.  082 E -02 

50 

1.  078E-02 

i.  078E-02 

70 

5.  550E-05 

5.  550E-05 

100 

1.  969E-08 

1.  969E-08 

by  the  values  of  N  in  Table  3  in  order  to  obtain  the  actual  size  distribution  at  any 
height  in  the  atmosphere. 

The  refractive  index  for  the  aerosol  particles  has  been  assumed  real  (that.  is. 
no  absorption)  for  Is  0.6  Mm.  For  *  >  o.  6  Mm  the  imaginary  part  of  the  refrac¬ 
tive  index  (that  is,  absorption)  is  assumed  to  increase  linearly  to  a  value  of  0. 1 
for  X  as  3  Mm.  This  model  for  the  refractive  index  is  based  on  measurements  by 
Vols  (19S7). 

The  total  numbers  (Table  3)  of  aerosol  particles  N  per  unit  volume  for  the 
"clear"  atmosphere  have  been  adjusted  so  that  the  total  extinction  coefficient  at 


wavelength  X  »  0.  55  nm  becomes  identical  to  the  values  used  in  the  atmospheric 
attenuation  model  by  Elterman  (1968,  1970)  at  each  altitude.  This  model  corres¬ 
ponds  to  an  atmospheric  condition  of  approximately  2  3  km  ground  visibility.  The 
"hazy"  model  is  identical  to  the  "clear"  model  for  altitudes  H  >  5  km.  For 
H  <  5  km,  the  total  number  of  aerosol  particles  increases  exponentially  to  reach 
a  value  at  ground  level  corresponding  to  a  ground  visibility  of  approximately  5  km 


3.  FUNDAMENTALS  OF  ATMOSPHERIC  TRANSMITTANCE 

Transmittance  of  radiation  in  the  atmosphere  is  complex  owing  to  the  de¬ 
pendence  of  scattering  and  absorption  coefficients  on  a  number  of  different  physical 
properties  of  the  atmosphere  (Goody,  1964).  In  general,  we  have  for  the  mono¬ 
chromatic  transmittance  of  radiation  along  a  path  in  the  atmosphere 

r  =  exp  (  -y  A  L )  (3) 

where  y  is  an  attenuation  coefficient  and  A  L  is  the  length  of  the  path  traversed  by 
the  radiation.  The  attenuation  coefficient  (y)  is  given  by 

y  =  a  +  k  (4) 

where  a  is  the  scattering  coefficient  and  k  is  the  absorption  coefficient.  Both  of 
these  coefficients  are  defined  in  km-1.  Equation  3  is  only  strictly  valid  when 
applied  to  monochromatic  radiation. 

The  scattering  and  absorption  coefficients  defined  in  Eq.  (4)  depend  on  both 
molecules  and  aerosols  in  the  atmospheric  path.  If  electromagnetic  radiation  is 
incident  on  a  molecule  or  an  aerosol,  a  portion  of  the  radiation  is  absorbed  and 
the  rest  is  scattered  in  all  directions.  Thus,  we  must  make  the  following  defini¬ 
tions  : 


a  *  a  +  o 

m  a 


(4a) 


k  +  k 

m  a 


(4b) 


where  the  subscripts,  m  and  a,  indicate  molecule  and  aerosol  respectively.  We 
will  endeavor  in  this  report  to  provide  quantitative  information  on  all  four  of  the 
quantities  defined  in  Eqs.  (4a)  and  (4b). 

The  molecular  scattering  coefficient  depends  only  on  the  number  density  of 
molecules  in  the  radiation  path,  whereas  the  molecular  absorption  coefficient  is 


i 


11 


a  function  of  not  only  the  amount  of  absorbing  gas,  but  also  the  local  temperature 
and  pressure  of  the  gas.  The  wavelength  dependence  of  molecular  (Rayleigh) 
scattering  is  very  nearly  am  ~  X"4  .  The  variation  of  the  molecular  absorption 
coefficient  with  wavelength  is  much  more  complicated,  being  a  highly  oscillatory 
function  of  wavelength  due  to  the  presence  of  numerous  molecular  absorption  band 
complexes.  These  band  complexes  result  for  the  most  part  from  minor  atmos¬ 
pheric  constituents.  The  responsible  molecules  are  (in  order  of  importance):  H20, 
COg,  O3,  N20,  CO,  02.  CH^  ,  N2.  All  of  these  molecules  except  H20  and  Og 
are  assumed  to  be  uniformly  mixed  by  volume. 

The  quantities  a&  and  ka  depend  on  the  number  density  and  size  distribution  of 
aerosols  as  well  as  on  their  complex  index  of  refraction.  The  amount  of  absorbed 
radiation  and  the  distribution  of  scattered  light  can  be  derived  theoretically  for 
spherical  particles  (and  some  other  simple  forms)  by  a  theory  developed  by  Mie 
(1908).  According  to  this  theory  the  distribution  of  radiation  with  scattering  angle 
^can  be  expressed  as: 


where  ij  and  i2  represent  the  two  polarized  components  of  the  scattered  radiation 
(van  deHulst,  1957).  They  are  functions  of  the  ratio  of  particle  size  to  wavelength, 
the  particle  refractive  index,  and  the  scattering  angle  <p.  I  is  the  intensity  of 
Incident  radiation.  If  there  are  N  particles  per  unit  scattering  volume  with  a  size 
distribution  n(r),  replacing  ij  and  ig  by  /  lj>2(r)n(r)  dr  will  yield  the  integrated 
angular  scattering  intensity  I(^)  from  the  unit  volume  of  the  particle  cloud  into  one 
steradian. 

e  4  w 

The  integral  /  I(^)d SI  gives  the  total  amount  of  scattered  radiation  which 

is  also  equal  to  I  •  . 

o  a 

Eq.  (5)  can  be  written  as: 

!(*)  ■  IQ  •  a  •  (6) 

where  4fy^)  is  the  normalized  scattering  function,  and  the  integral 

/4Vw  dfl  becomes  1. 
o 

Figure  26  shows  ^»(^)  for  molecular  scattering  (wav-length  independent)  and 
aerosol  scattering  for  two  wavelengths  for  the  aeroBol  model  described  in 
Section  2.3.  An  absorption  coefficient,  kft,  for  aerosol  particles  can  be 
defined  in  the  same  manner  as  the  scattering  coefficient,^.  The  absorption  and 
scattering  coefficients  due  to  aerosol  distributions  are  usually  smooth  functions 
of  wavelength. 
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The  molecular  absorption  coefficient  for  a  single  spectral  line  can  be  given  to 
good  accuracy  in  the  lower  50  km  of  the  atmosphere  by 


k  (v> 
m 


*  ft"  ‘  "0)2  +  Q'2] 


where  S  is  the  line  intensity,  a  is  the  line  half-width,  vQ  is  the  frequency  of  the 
line  center  and  v  .s  the  frequency  at  which  the  absorption  coefficient  is  required. 
See  the  Appendix  for  the  relationship  between  v  and  X.  The  half-width  depends  on 
the  pressure  and  temperature  approximately  as  follows: 


P 

ao  F- 
o 


The  line  intensity  also  depends  on  the  temperature  through  the  Boltzmann  factor 
and  the  partition  function  (Goody,  1964  and  Penner,  1959).  Above  50  km,  Doppler 
broadening  (Goody,  1964)  should  be  considered. 

Let  us  consider  a  plane  parallel  atmosphere  and  a  frequency  at  which  several 
lines  belonging  to  one  or  more  molecular  constituents  contribute  to  the  absorption 
coefficient.  The  molecular  absorption  coefficient  in  any  one  of  the  (j)  layers  will 
be  given  by: 


kmj  ‘  I  kmji 


where  the  (1)  summation  is  over  all  molecular  absorption  lines  due  to  all  absorbing 
species  which  are  close  enough  to  the  frequency  v  to  contribute  significantly  to  the 
total  molecular  absorption  coefficient  defined  In  Eq.  (9).  Introducing  the  molecular 
scattering  ( am ),  and  the  absorption  (fc^)  and  scattering  (oa)  due  to  aerosols  as 
defined  in  Eqs.  (4a)  and  (4b),  we  have  for  the  total  attenuation  coefficient  in  the 
(jth)  layer: 

Vkmj  +  CTmJ+llaj  +  aar  (10) 


Thus,  according  to  Eq,  (3),  the  transmittance  of  an  atmosphere  having  (J)  layers 
is  given  by: 

r(v)  •  exp  ALj)*  <,]l 
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It  should  be  noted  from  Eq.  (11)  that  AL  does  not  necessarily  denote  vertical  dis¬ 
tance.  The  problem  of  slant  path  transmittance  is  easily  solved  by  writing 
ALj  ■  (&Zj)  sec  6  where  represents  the  increment  of  vertical  distance.  In 
this  case,  Eq.  (11)  becomes 

tm(v)  -  [r(t/)]M  (12) 

where  M  =  sec  0  is  the  secant  of  the  angle  between  the  vertical  and  the  direction  of 
observation.  Thus  it  becomes  very  easy  to  convert  vertical  transmittance  values 
to  slant  path  transmittance.  As  0  approaches  90°,  the  result  becomes  meaningless 
due  to  the  plane  parallel  atmosphere  assumption  (see  Section  7). 

All  of  the  above  discussion  is  valid  only  at  a  single  wavelength  (or  very  nar¬ 
row  wavelength  interval  which  is  known  to  be  small  compared  with  the  width  of  any 
absorption  lines  in  the'vicinity).  A  sufficiently  narrow  spectral  Interval  for  most 
purposes  is  0.001  wavenumber  (cm~*).  The  process  which  so  limits  the  validity 
of  the  above  relations  is  molecular  absorption.  As  noted  above,  k^fv)  is  a  highly 
oscillatory  function  of  frequency  whereas  ka(v),  and  <r  (v)  are  much  more 

slowly  varying  functions. 

When  it  is  desired  to  calculate  the  atmospheric  transmittance  over  a  large 
spectral  interval,  it  is  necessary  to  average  Eq.  (11)  or  Eq.  (12)  over  the  desired 
spectral  interval  (see  Eq.  (14)].  When  many  absorption  lines  are  present,  this 
can  be  a  particularly  tedious  job  (even  on  a  computer).  Various  techniques  have 
been  developed  to  treat  this  problem  and  are  discussed  in  Section  5. 

4.  ATMOSPHERIC  TRANSMITTANCE  AT  SELECTED  FREQUENCIES 

4.1  List  of  Laser  Frequencies 

It  is  necessary  to  limit  the  number  of  specific  frequencies  at  which  we  quote 
atmospheric  absorption  and  scattering  coefficients.  We  have  therefore  chosen  the 
list  of  specific  laser  wavelengths  (Table  4)  at  which  values  of  these  coefficients 
are  given  in  Table  5. 

4.2  Scattering  and  Absorption  Coefficients 

Each  atmospheric  model  defined  in  Section  2  has  been  used  as  a  basis  for 
calculating  the  absorption  and  scattering  coefficients  given  in  Table  5  for  each 
laser  wavelength.  Numerical  values  for  the  coefficients  were  computed  for  the  0, 

1,  2,  etc. ,  km  levels  and  then  logarithmically  interpolated  for  each  atmospheric 
layer.  The  aerosol  attenuation  coefficients  given  In  Table  5,  which  correspond  to 
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Table  4.  Laser  Wavelengths^ 


0. 3371  urn 

Nitrogen 

0. 4880  /urn 

Argon 

0.  5145  ^m 

Argon 

0.  6328  gm 

Helium -Neon 

0.  6943  jum 

Ruby 

0. 86  /urn 

Gallium  Arsenide 

1. 06  jum 

Neodymium  in  Glass 

1.  536  jum 

Erbium  in  Glass 

3.  39225/um 

Helium  Neon 

10.591  Mm 

Carbon  Dioxide 

27.  9  pm 

Water  Vapor 

337  utm 

Hydrogen  Cyanide 

the  clear  and  hazy  atmospheres  defined  in  Section  2.  3,  are  Identical  for  all  geo  ¬ 
graphical  and  seasonal  models.  In  Table  5,  the  four  extinction  coefficients  are 
listed  as  a  function  of  altitude  for  each  model:  ^  and  <7m  are  the  molecular  ab¬ 
sorption  and  scattering  coefficients,  and  and  oa  are  the  aerosol  absorption  and 
scattering  coefficients.  If  the  total  attenuation  coefficient  is  required,  the  four 
coefficients  can  be  summed: 


+  k„  +  a 

a  a 


(13) 


4.3  Interpolation  Procedures 

In  general,  it  is  permissible  to  interpolate  linearly  with  respect  to  wavelength 
between  tabulated  values  in  Table  5  except  for  molecular  absorption  coefficient 
(k^).  This  is  true  because  of  the  smooth  variation  of  scattering  coefficients  and 
of  aerosol  absorption  coefficients  with  wavelength,  and  the  rapid  spectral  variations 
in  molecular  absorption  coefficients.  The  tables  have  been  constructed  so  that 
linear  Interpolation  with  respect  to  height  is  expected  to  give  correct  values  to 
±8  percent.  If  a  given  atmospheric  condition  differs  from  one  of  the  ten  models  of 
Section  2,  larger  errors  may  result. 

Again  it  should  be  remembered  that  one  should  never  attempt  to  interpolate 
with  respect  to  wavelength  between  molecular  absorption  coefficient  values  0^ > 
given  in  Table  5, 

*  These  lasers  may  emit  over  a  range  of  wavelengths.  However,  the  computed 
values  in  Table  5  were  derived  for  the  specific  wavelength  values  of  Table  4. 


ile  5.  Values  of  Attenuation  Coefficient/km  as  a  Function  of  Altitude  for  Each  Laser  Wavelength  and  Each 
lospheric  Model  in  Section  2  (kjj;  =  molecular  absorption,  cr^  =  molecular  scattering.  =  aerosol 
orption,  o_  -  aerosol  scattering) 
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Having  obtained  the  attenuation  coefficient,  7  [Eq.  (13)],  the  laser  transmit¬ 
tance  for  a  horizontal  path  of  length,  L,  at  height,  h,  is  given  by  r  =  exp  (-7jL) 
when  the  suffix  j,  refers  to  the  layer  at  height,  h.  For  slant  paths  between  two 
altitudes,  Eq.  (11)  should  be  used  with  the  appropriate  correction  for  zenith  angle. 


5.  ATMOSPHERIC  TRANSMITTANCE  AT  LOW  RESOLUTION  FROM 
0.2S  TO  25  n  m 


The  calculations  presented  in  Section  4  were  made  at  infinite  resolution 
(monochromatic  calculations).  In  practice  it  is  impossible  to  measure  radiation 
transmitted  at  a  single  frequency.  Instead  one  measures  the  transmittance  (v) 
averaged  over  the  spectral  interval,  Av  ,  accepted  by  the  receiver,  as  indicated 
in  the  equation 


V 


**>  m~Sv 


hf 


T  ( V )  dv. 


(14) 


where  v  is  the  central  frequency  in  the  interval,  Ar.  Consequently  for  many  appli¬ 
cations  one  is  interested  in  knowing  the  transmittance  of  the  atmosphere  averaged 
over  a  relatively  wide  spectral  interval,  that  is,  for  low  resolution. 

Thus  the  term  transmittance  is  somewhat  ambiguous  unless  it  is  qualified  by 
some  indication  of  the  spectral  resolution,  Ar,  over  which  it  is  averaged.'  This  is 
particularly  true  in  the  case  of  molecular  absorption,  oince  the  absorption  coeffi¬ 
cient  kj^  is  a  rapidly  varying  function  of  frequency.  It  is  because  of  the  rapid 
variation  of  with  frequency  that  the  averaged  transmittance  r  does  not,  in 
general,  obey  the  simple  exponential  law  given  in  Eq.  (3).  That  is, 


TAi'(i/>  "  ST  /-P  d,/  *  exP  (•K(«/)SL)  (15) 

where  km  represents  the  net  monochromatic  molecular  absorption  coefficient  as 
given  in  Eq.  (9)  and  where  K  is  an  average  absorption  coefficient,  which  cannot  in 
most  cases  be  defined  when  Ar  is  much  greater  than  the  half  width  of  a  spectral 
line  (see  Eq.  (8)] . 

On  the  other  hand  the  molecular  scattering  coefficient  (om)  and  the  aerosol 
scattering  and  absorption  coefficients  (oa  and  kg,  see  Eqs.  (4a)  and  (4b)]  are  slow¬ 
ly  varying  functions  of  frequency,  and  the  average  transmittance  (which  can  easily 
be  obtained  by  Interpolation  from  the  monochromatic  values,  see  Figures  21  and  22) 
obeys  the  simple  exponential  law  (Eq.  (3)]  provided  only  the  direct  transmitted 
beam  is  being  observed  (see  Zuev  et  al,  1967). 
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There  are  four  basic  approaches  to  obtaining  low  resolution  transmittance 
values  for  a  given  path  through  the  atmosphere  due  to  molecular  absorption.  These 
are:  (1)  direct  measurements  over  the  required  path,  (2)  measurements  in  the 
laboratory  .under  simulated  conditions,  (3)  line-by-line  (monochromatic)  calcula¬ 
tions  based  on  detailed  knowledge  of  spectroscopic  line  parameters  (see  Section  3) 
which  are  then  averaged  over  the  required  spectral  interval,  and  (4)  calculations 
based  on  band  model  techniques  (which  use  available  laboratory  and/or  field  trans¬ 
mittance  measurements  or  actual  line  data  as  a  basis). 

From  the  point  of  view  of  computations,  method  3  involves  a  considerable 
amount  of  work  and  computer  time,  and  consequently  method  4  has  been  used  most 
frequently. 

In  this  section  we  will  briefly  discuss  the  problem  of  continuum  absorption  in 
the  atmospheric  window  region  (Section  5.  1).  Then  we  will  summarize  the  appli¬ 
cation  of  band  models  to  molecular  absorption  calculations  (section  5.2)  and 
present  a  detailed  description  of  a  graphical  prediction  scheme  for  determining 
the  transmittance  of  the  atmosphere  due  to  line  absorption,  continuum  absorption, 
and  scattering  (average  over  20  cm*1)  for  a  given  atmospheric  path.  First,  the 
reader  is  provided  with  figures  for  determining  the  relevant  absorber  amounts  in 
the  required  atmospheric  path  (Section  5.  3),  which  are  then  used  in  conjunction 
with  the  transmittance  charts  (Section  5.4)  to  determine  (separately)  the  transmit¬ 
tance  due  to;  (1)  molecular  (line)  absorption  (H20,  03,  and  the  uniformly  mixed 
bases),  (2)  molecular  (continuum)  absorption  (H20,  N2),  (3)  molecular  (Rayleigh) 
scattering,  and  (4)  aerosol  extinction  (absorption  and  scattering),  as  a  function  of 
frequency  (wavelength). 

Finally  die  reader  has  to  multiply  together  the  individual  transmittance s  at 
each  wavelength  in  order  to  determine  the  total  transmittance  at  that  wavelength. 
That  is: 

T^fv)  (Total)  ■  (p)  (line  absorptton)  x  (v)  (continuum  absorption) 

y  “^(v)  (Rayleigh)  x  t^v)  (aerosol)  . 

5.1  Continuum  Absorption  In  the  Atmoapherin  f  tn<io>v  Region* 

There  are  a  number  of  spectral  regions  of  weak  or  negligible  absorption  by 
atmospheric  molecules  in  the  infrared  which  are  generally  called  atmospheric 
windows.  Various  attempts  have  been  made  to  determine  the  transparency  of  these 
window  regions  (Yates  and  Taylor,  1960;  and  Streets,  1968),  Yates  and  Taylor 
(1960)  concluded  that  the  opacity  in  the  near  infrared  windows  is  limited  at  sea 
level  by  aerosol  scattering.  This  conclusion  is  consistent  with  the  results  which 
are  derived  from  our  aerosol  chart  presented  in  Section  5. 4 
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As  there  are  no  laboratory  data  presently  available  from  which  the  molecular 
absorption  between  absorption  bands  can  be  inferred  (with  two  important  excep¬ 
tions),  we  will  here  assume  that  the  gaps  in  the  charts  presented  in  Figures  16 
through  18  are  real  gaps  of  zero  molecular  absorption.  Thus,  the  attenuation  in 
these  regions  is  dependent  entirely  on  molecular  and  aerosol  scattering. 

However,  laboratory  absorption  data  are  available  for  the  water  vapor  con¬ 
tinuum  (Burch,  1970;  and  McCoy  et  al,  1969),  from  8- 13 Mm  and  also  for  the  pres¬ 
sure  induced  (continuous  absorption)  band  of  nitrogen  (Bosomworth  and  Gush,  1965) 
centered  near  4.  3  pm.  The  strongest  part  of  this  nitrogen  band  is  not  too  import¬ 
ant  as  it  is  centered  beneath  the  very  strong  4.  3  /um  C02  band.  However,  the 
nitrogen  absorption  extends  out  beyond  the  region  of  C02  absorption  near  4.  0  pm 
and  becomes  predominant,  accounting  for  20  percent  absorption  in  a  10  km  path 
at  sea  level. 

Due  to  the  continuous  nature  of  these  absorbing  regions,  the  transmittance 
follows  a  simple  exponential  law  which  is  represented  by  the  scale  in  Figures  19 
and  1’0. 

5.2  Outline  of  Hand  Model  Technique.* 

Basically,  a  band  model  assumes  an  array  of  lines  having  chosen  intensities, 
half-widths  and  epacings,  which  can  be  adjusted  to  represent  the  line  structure  in 
some  part  of  a  reaJ  band.  For  a  particular  band  model  the  mean  transmittance  can 
be  represented  by  a  mathematical  expression  (transmittance  function);  for  example 

rAl/(r)  »/ (C(r),  AL,P)  (16) 

expressed  in  terms  of  pressure,  P,  effective  path  length  (or  absorber  concentra¬ 
tion)  £L,  and  one  or  more  frequency  dependent  absorption  coefficients  (C(i')'s). 

Several  band  models  have  been  developed  (see  Goody,  1952  and  1964;  Flsasser, 
1938;  Plass,  1958;  King,  1959  and  1964;  and  Wyatt  et  al,  1964);  the  Flsasser  mode) 
(1942)  and  the  Goody  model  (1952  and  1964)  being  most  well  known.  The  Flsasser 
band  model  consists  of  an  array  of  equally  spaced  identical  lines  of  Forentz  shape 
(see  Fq,  (7) J .  This  model  has  been  applied  to  absorption  bands  which  have  a 
regular  line  structure,  for  example,  to  some  bands  of  C02,  N2O,  CO,  CH4  and 
02<  The  Goody  model,  on  the  other  hand,  assumes  that  the  band  is  composed  of 
spectral  lines  with  an  exponential  intensity  distribution  and  with  random  spacing 
between  lines.  Again  the  Mnes  are  assumed  to  have  a  Forentz  shape.  This  model 
has  generally  been  applied  to  bands  which  have  an  irregular  line  structure,  for 
example,  to  some  H^O  and  bands. 

In  practice  the  wavelength  dependent  coefficients  are  determined  for  each 
absorbing  gas  separately  using  laboratory  transmittance  values  measured  under 
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known  condition*,  that  ia  by  solving  Eq.  (16)  for  the  C(v)'s  at  each  v  for  known 
values  of  t(v),  AL  and  P, 

The  absorption  coefficient*  obtained  in  this  way  are  then  used  in  the  band 
model  transmittance  function  (Eq.  (16))  to  determine  r( v)  as  a  function  of  frequency 
v  for  other  value*  of  AL  and  P  for  each  absorber.  Finally,  the  total  mean  trans¬ 
mittance  for  molecular  absorption  is  given  by  the  product  of  the  mean  transmit¬ 
tance*  of  the  individual  absorbers  at  each  frequency. 

V"  ( v)  (molecular  absorption)  ■  T^y(i/)  (H^O)  *  ^^(r)  (C02) 

*  ?Av(N20)  *  <03>  etCl  (17) 


Exact  analytical  expressions  have  been  obtained  for  most  of  the  band  models. 
However,  they  are  sometimes  difficult  to  use  and  simpler  approximations  have 
been  found  to  apply  in  two  limiting  conditions,  which  are  common  to  all  band 
models.  These  simpler  expressions  are  the  well  known  "weak  line"  and  "strong 
line"  approximations  (sec  Goody,  1964;  and  Plass,  1958)  for  which  the  transmit¬ 
tance  is  a  function  of  the  absorber  amount,  AL,  and  the  product  of  the  pressure  and 
absorber  amount,  PAL,  respectively  (fora  given  temperature). 

The  weak  line  approximation,  which  corresponds  to  the  exponential  law,  is 
valid  when  the  absorption  is  small  at  the  line  centers  (generally  for  high  pressures 
and  low  absorber  amounts).  Unfortunately  this  case  is  rarely  applicable  to  condi¬ 
tions  existing  in  the  terrestrial  atmosphere. 

The  strong  line  approximation  is  applicable  where  the  lines  are  completely 
absorbing  at  their  centers;  the  effect  of  increasing  the  amount  of  absorber  is  then 
confined  to  the  edges  or  wings  of  the  lines. 

The  regions  of  validity  of  the  strong  and  weak  line  approximations  for  the 
Elsasser  and  Goody  models  are  discussed  in  Plass  (1958), 

For  practical  purposes  most  problems  fall  in  either  the  strong  line  approxima¬ 
tion  region  or  the  Intermediate  regime,  it  is  perhaps  interesting  to  compare  the 
functional  form  of  the  strong  line  approximation  for  some  of  the  models. 

For  the  Elsasser  Model 


tAv<1')  e  1  '  erf  (CM  ^LPJ 


where  erf  la  the  error  function 


1/2 


erf  (z) 


=  _2_  f 

Jn~  J  z 


-l 

e  1  dt 


(18) 


and  C(v)  is  an  average 


\Jn~  ^  z 

absorption  coefficient  which  can  be  written  in  terms  of  average  values  of  the  mole¬ 
cular  line  constants  defined  In  Section  3. 


The  strong  line  approximation  for  the  Goody  model  is  given  by 


T6vM  "  exP  {[C'^)  ALP]1/2}  (19) 

where  C'(r)  is  defined  in  a  similar  manner  to  C(v)  above. 

King  (1959)  attempted  to  generalize  the  above  expressions  by  writing 

tAi/(v)  =/ (c"(i/)  ALPnj  (20) 

so  that  when  n  is  set  equal  to  zero  or  unity,  the  average  transmittance  agrees 
with  the  weak  or  strong  line  approximations  respectively. 

In  this  report  we  have  determined  from  laboratory  and  synthetic  transmittance 
data,*  three  empirical  transmittance  functions  based  on  Eq.  (20)  for  (1)  H2O, 

(2)  O3,  and  (3)  the  combined  contributions  of  the  uniformly  mixed  gases.  It  was 
found  that  Eq.  (2  0)  gave  the  best  fit  to  laboratory  and  theoretical  data  over  a  wide 
spectral  interval  when:  (1)  n  =  0.  9  for  HgO,  (2)  n  =  0.  4  for  O3,  and  (3)  n  =  0.  75 
for  the  uniformly  mixed  gases. 

The  corresponding  empirical  transmittance  functions  were  found  to  give  better 
agreement  with  laboratory  and  synthetic  transmittance  data  than  the  commonly 
used  band  models  over  a  wide  range  of  pressures  and  absorber  amounts. 

The  procedure  used  for  determining  the  parameter  n  and  function,  J,  defined 
in  Eq.  (20)  will  be  briefly  outlined  below.  By  taking  the  logarithm  of  the  inverse 
of  Eq.  (20),  it  wi.1  be  seen  that 

n  log  P  +  log  AL  -  log/"1^^)]  -  log  C"  (v>.  (21) 

Thus  for  a  given  frequency  v  and  fixed  values  of  average  transmittance 
t  (i/),  the  righthand  side  of  Eq.  (20)  becomes  a  constant. 

A  mean  value  for  n  was  determined  from  Eq.  (21)  for  a  wide  range  of  fre¬ 
quencies  and  several  values  of  T.  Then  for  each  frequency,  r  was  plotted  against 
log  AL  Pn,  and  the  curves  superimposed  and  the  best  mean  curve  determined. 

The  "mean"  curve  thus  obtained  constitutes  the  empirical  transmittance  function 
and  is  displayed  as  the  transmittance  scale  and  associated  scaling  factor, 
log  (ALFn),  in  the  prediction  charts  given  in  Figures  16  through  18. 

So  far  the  discussion  has  been  restricted  to  constant  pressure  paths.  In  the 
following  section  we  will  discuss  the  application  of  band  models  to  variable  pressure 
paths. 

♦Synthetic  spectra  here  refers  to  monochromatic  (line  by  line)  transmittance 
calculations  degraded  in  resolution  to  20  cm*’1. 


/  -  V 


32 


S.3  Equivalent  Constant  Pressure  Paths  and  Absorber  Amounts 

Atmospheric  paths  may  be  divided  into  two  types,  namely,  constant  pressure 
paths  (that  is,  horizontal  paths  of  up  to  about  100  km)  and  slant  or  variable  pres¬ 
sure  paths.  The  extent  of  the  former  types  of  path  is  limited  due  to  the  curvature 
of  the  earth.  Since  the  transmittance  functions-  for  constant  pressure  paths  are 
simpler  to  handle  than  those  for  variable  pressure  paths,  the  problem  is  simplified 
by  reducing  the  slant  path  to  an  "equivalent  constant  pressure  path";  that  is,  to  a 
constant  pressure  path  which,  ideally,  has  the  same  transmittance  at  each  wave¬ 
length  as  does  the  actual  slant  path. 

In  this  section  we  will  discuss  the  construction  and  application  of  charts  for 
determining  the  equivalent  sea  level  path  quantities  for  molecular  absorption  (in¬ 
cluding  continuum  absorption  in  the  4  pm  and  8  -13  pm  region),  molecular  scatter  ¬ 
ing  and  aerosol  extinction  based  on  the  model  atmospheres  presented  in  Section  2. 

5.3.1  MOLECULAR  ABSORPTION 


It  has  been  assumed  that  for  most  practical  purposes,  the  average  trans¬ 
mittance  is  dependent  on  the  product  ALPn  where  AL  is  the  absorber  amount  and 
P  is  the  pressure. 

Here  we  will  define  the  equivalent  sea  level  amounts  for  each  of  the  five  model 
atmospheres  given  in  Section  2  as  the  value  of  the  above  product  (reduced  to  3TP) 
averaged  over  a  given  atmospheric  path  for  a  particular  absorber,  that  is 

A  Lq  where  tne  subscript  (o)  refers  to  STP. 


Since  II20  and  03  are  distributed  nonuniformly  in  the  atmosphere,  we  must 
consider  these  gases  separately  from  the  uniformly  mixed  gases. 

For  water  vapor,  n  is  set  equal  to  0.  9  and  AL0  (^rj  a  Rw  9 


where 


w(z)  is  the  water  vapor  concentration  at  altitude,  z,  in  units  of  gm  cm'2/km 
(obtained  from  Table  1),  and  K  is  the  range  in  km.  Since  w(z)  is  given  in  mass 
units  there  is  no  temperature  or  pressure  correction  necessary  to  reduce  this 


quantity  to  STP. 


The  variation  with  altitude  of  w(z)  and J“  w(z)^pr^ 


dz 


(denoted  by  wh  and  wy  respectively)  is  shown  in  Figures  2  and  3.  The  subscripts, 
h  and  v,  refer  to  horizontal  and  vertical  paths  respectively. 

For  the  uniformly  mixed  gases,  n  is  set  equal  to  0.75,  but  the  absorber  con¬ 
centration,  AL0  (reduced  to  STP)  for  a  horizontal  path  of  range  R  at  altitude  z 

is  given  by  AL0  ■  cR  &)(■*)  where  c  is  the  fractional  concentration  by  volume  of 

the  absorbing  gas  (see  Table  2  and  Appendix).  The  terms  P  and  T  refer  to  the 
pressure  and  tempera  tore  at  altitude  z.  Therefore,  we  have  in  this  case 
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^Lo(?0  (-^)  •  The  variation  with  altitude  of  (~^j'  ?5  ^  and 

/•°°/pvl.75  /T  \ 

J  \P0/  XT’/  7'  denoted  by  and  fy  respectively),  is  shown  in  Figures 


4  and  5. 


vO.  4 


For  ozone,  n  is  set  equal  to  0.  4  (see  Section  5.  2),  andALQ^-J  «  u(z) 
where  u(z)  is  the  ozone  concentration  in  units  of  (cm  -  atm)/km  at  STP  (obtained 

from  Table  1).  The  variation  with  altitude  of  u(z)^-J‘  4  and  J  u(z)^-J  4dz 
(denoted  by  uh  and  uy  respectively)  is  shown  in  Figures  6  and  7. 

5.  3.2  NITROGFN  CONTINUUM  (4  um  REGION) 

The  absorption  due  to  the  nitrogen  collision  induced  band  is  proportional  to 
the  square  of  the  pressure.  The  effective  absorber  amount  (reduced  to  STP)  for 

a  horizontal  path  of  range  R  at  altitude  z  is  cR 


(£)  (t)  •  The 

rm) 


altitude  off^-^  l-=?Jfor  horizontal  paths  and 
(denoted  by  n^  and  ny  respectively)  is  shown  in  Figures  8  and  0 
5.  3.  3  WATER  VAPOR  CONTINUUM 


variation  with 

dz  for  vertical  paths 


Laboratory  measurements  of  the  water  vapor  continuum  indicate  that  self¬ 
broadening  due  to  water  vapor  molecules  is  much  more  important  than  broadening 
due  to  other  atmospheric  molecules  (Burch,  1970).  The  appropriate  equivalent 
sea  level  absorber  amount  for  the  water  vapor  continuum  is 


bh  “  ^  [Ph2°  +  °*  005  (P  '  PH2o)]  w(z) 

for  horizontal  paths,  where  P  is  the  total  pressure  and  Ph20  is  the  partial  pressure 
of  water  vapor  (Ph20  1  4.56  xl0‘6  w(z)T(z)  atm,  where  the  constant  is  given  by  the 
ratio  of  the  universal  gas  constant  to  the  molecular  weight  of  water  vapor)  and  w(z) 
is  the  number  of  (gm  cm"2/km)  of  water  vapor  at  altitude  z  obtained  from  Table  1. 
The  equivalent  sea  level  amount  for  vertical  paths  is  given  by 

w(z)dz  . 
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The  variation  with  altitude  of  these  functions,  denoted  by  and  by  respectively, 
is  shown  in  Figures  10  and  11. 

5.3.4  MOLECULAR  SCATTERING 


Molecular  scattering  is  a  function  of  atmospheric  density  p(z)  (Section  3). 

For  each  of  the  five  model  atmospheres  given  in  Section  2,  Figures  12  and  13  show 

r 

the  variation  with  altitude  z  of  p(z)/p  and  /  p(z)/pQ  dz  (denoted  by  and  m 

Jz  v'  **  v 

respectively)  where  p  is  the  density  of  air  at  STP. 


5.  3.  5  AEROSOL  EXTINCTION 


Aerosol  extinction  is  a  funci  i  of  the  height  variation  of  the  aerosol  number 
density  N(z)  for  a  given  size  distribution  function  (Section  2.  3).  Figures  14  and  15 


N(z) 

show  the  variation  with  altitude  of  -rj-  and 

Nc  -z 


r 


N(z)dz 


~ - (denoted  by  a,  and  a 

iNq  n  v 

respectively)  for  the  two  hazy  models  given  in  Section  2.  3  where  N0  is  the  number 

density  (cm*3)  of  aerosols  at  sea  level  for  a  visibility  of  23  km. 


5.3.6  UTILIZATION  OF  EQUIVALENT  PATH  FIGURES 


For  a  given  atmospheric  path,  let  us  consider  how  to  obtain  the  equivalent 
sea  level  quantities  for  HgO.Og,  the  uniformly  mixed  gases,  the  nitrogen  and 
water  continua,  molecular  scattering  and  aerosol  extinction  which  will  be  denoted 
by  W,  U,  L,  N,  B,  M,  and  A  respectively.  At  this  point  the  reader  need  not  be 
concerned  with  the  units  of  these  quantities  which  were  chosen  to  facilitate  the  use 
of  the  prediction  charts  given  in  Section  5.4.  The  subject  of  units  will  be  discussed 
separately  in  the  Appendix. 

Figures  2  through  15  show  the  variation  with  altitude  of  the  relevant  quantities 
for  both  horizontal  and  vertical  paths  (distinguished  by  suffixes  h  and  v  respec¬ 
tively).  The  use  of  these  figures  can  best  be  illustrated  by  the  following  examples: 

5.  3.  6.  1  Horizontal  Paths 


To  determine  the  equivalent  sea  level  quantities  for  a  horizontal  path  of 
length  H  (<  100  km)  at  an  altitude  z,  simply  multiply  each  of  the  quantities  w^, 
uh’  nh*  mh  and  ah  Stained  from  Figures  2,  4,  6,  8,  10,  12,  14  for  the  ap¬ 
propriate  altitude  z),  by  R  (in  tons).  For  example  in  the  case  of  molecular  line 
absorption  by  water  vapor,  W  is  given  by 


W  -  wh(z)R. 


(21) 
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5.  3. 6. 2  Slant  Paths 

(1)  The  equivalent  sea  level  quantities  (W,  L»  U,  N,  B,  M  and  A)  for  a  slant 
path  from  altitude  z  to  the  top  of  the  atmosphere  at  zenith  angle  8  are  determined 
by  multiplying  each  of  the  quantities  wy,  iy,  uy,  ny,  by,  my,  and  ay  (obtained 
from  Figures  3,  5,  7,  9,  11,  13,  15  for  the  appropriate  altitude  z)  by  sec  8  (where 
8  <  80°).  For  example,  in  the  case  of  molecular  line  absorption  by  water  vapor 

W  «  wy(z)  sec  9  .  (22) 

For  8  >  80°,  sec  8  should  be  replaced  by  M(0)  given  in  Figure  32  (see Section  7). 

(2)  For  a  slant  path  of  range  R  between  altitudes  z  j  and  z2>  the  value  of  W  for 
molecular  line  absorption  by  water  vapor  is  given  by 

W  •  [wy  (Zl)  -  Wy  <z2>]  X  Tpriq  *  <23> 

Similar  expressions  hold  for  L,  U,  N,  B,  M  and  A. 

Corrections  can  also  be  applied  to  account  for  atmospheric  refraction  and 
earth  curvature  and  these  will  be  briefly  summarized  in  Section  7. 

Having  now  determined  W  (gm.  cm”2),  U  (cm-atm  at  STP),  L  (km),  N  (km), 

B  (gm  cm-2),  M  (km)  and  A  (km)  one  can  use  the  prediction  charts  given  in 
Figures  16  through  22  (Section  5,4). 

5,4  Description  and  Use  of  Prediction  Charts 

The  purpose  of  this  section  is  to  present  a  simple  method  of  predicting  atmos¬ 
pheric  transmittance  (at  low  resolution)  which  is  applicable  over  a  wide  spectral 
interval  and  for  a  wide  range  of  atmospheric  paths.  Of  the  methods  available,  the 
one  originally  devised  by  Altshuler  (1961)  is  the  most  comprehensive  that  meets 
the  above  requirements,  although  based  on  fairly  old  experimental  data.  In  view 
of  this,  an  attempt  has  been  made  here  to  further  simplify  Altshuler's  method  and 
to  update  it  by  (1)  using  more  recent  laboratory  and  lleld  transmittance  data  (Burch 
et  al,  1962,  1964,  1965a,  1965b,  1965c,  1965d,  1967,  and  1968)  and  (2)  making 
some  line  by  line  calculations  for  molecular  absorption  (as  discussed  in  Section 
5.  4)  and  averaging  over  the  spectral  interval  (bv  ■  20  cm"*). 

Charts  for  calculating  the  transmittance  due  to  molecular  line  absorption  by 
water  vapor  (HgO),  ozone  (Og),  the  uniformly  mixed  gases  (COg.  NgO,  CH4,  CO, 
02),  the  nitrogen  4  gm  continuum,  the  water  vapor  8-13  gm  continuum,  molecular 
scattering,  and  aerosol  extinction  and  absorption  are  presented  in  Figures  16 
through  22.  It  will  be  noted  that  the  effects  of  the  different  uniformly  mixed  gases 
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have  been  combined  into  one  chart  (Figure  1?)  (assuming  the  mixing  ratios  given 
in  Table  2). 

A  description  and  outline  of  the  use  of  the  charts  will  now  be  given. 

Each  figure  shows  the  variation  of  an  attenuation  coefficient  with  respect  to 
wavelength,  relative  to  a  horizontal  datum  line  for  a  given  effective  sea  level  path 
length  or  absorber  amount.  Also  shown  in  each  figure  is  a  vertical  transmittance 
scale  associated  with  a  set  of  scaling  factors  (for  the  absorber  amount).  In  order 
to  use  these  charts  efficiently,  it  is  necessary  to  trace  the  transmittance  scale  and 
associated  scaling  factors  for  each  of  the  figures  on  to  some  transparent  paper. 

If  the  transmittance  scale  is  now  placed  on  its  appropriate  chart  so  that  the 
scaling  factor  1  rests  on  the  datum  line  (as  shown  in  Figure  16(b)),  then  one  is  in  a 
position  to  read  off  the  transmittance  (for  1  gm.  cm'2  of  H20)  as  a  function  of 
wavelength.  Simply  move  the  chart  horizontally  until  the  transmittance  scale  is 
vertically  above  the  required  wavelength  and  read  off  the  transmittance  where  the 
scale  crosses  the  curve.  For  example,  in  Figure  16(b)  tne  scale  is  set  at  21 
microns  and  indicates  a  transmittance  of  0.  10  for  W  «  1  gm  cm'2  of  H20.  *  If  the 
transmittance  is  required  for  a  different  amount  of  absorber,  for  example, 

W  *  10  gm  cm*2,  just  displace  the  transmittance  scale  vertically  until  the  scaling 
factor  10  is  coincident  with  the  datum  line.  Then  move  the  scale  horizontally  to 
the  required  wavelength  and  read  off  the  transmittance  as  described  above. 

The  same  procedure  applies  to  Figures  17  through  22,  25  and  33,  where  it 
will  be  noted  the  datum  line  for  scaling  the  absorber  amounts  is  expressed  as 
1  (cm-atm)gTp  for  Og  (Figures  18  and  25),  1  gm  cm'2  for  the  water  vapor  con¬ 
tinuum  (Figure  20),  and  1  km  for  each  of  the  other  charts. 

Thus  by  displacing  the  transmittance  scale  horizontally  according  to  wave¬ 
length  and  vertically  according  to  the  appropriate  equivalent  sea  level  quantity  (W, 

U,  L,  N,  B,  M  or  A)  for  Figures  16  through  22  respectively,  one  can  construct  a 
low  resolution  transmittance  spectrum  (for  &v  >20  cm'1). 

A  method  for  calculating  the  equivalent  sea  level  absorber  amounts  W,  U,  L, 

N,  B,  M  and  A  respectively  (in  the  appropriate  units)  for  a  given  atmospheric  path 
has  been  described  in  Section  5.  3  based  on  the  ten  model  atmospheres  defined  in 
Section  2.  In  the  case  where  the  relevant  meteorological  parameters  are  known 
for  the  atmospheric  path,  these  quantities  can  be  calculated  by  an  alternative 
method  described  in  the  Appendix. 

5.3  Abaorption  in  the  Region  0. 25  to  0.75  itm 

The  only  atmospheric  molecule  considered  in  this  report  in  the  spectral  region 
from  0.25  to  0.70pm  is  ozone.  Weak  absorption  bands  due  to  water  vapor  and  oxygen 

*This  corresponds  to  the  amount  of  water  vapor  in  a  1  km  path  at  sea  level  for 
T  ■  16°C  and  78  percent  relative  humidity  (see  Appendix). 
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also  appear  in  this  region  but  because  of  a  lack  of  reliable  data  have  not  been  in¬ 
cluded  here.  Since  the  absorption  coefficient  due  to  ozone  in  this  wavelength  range 
is  independent  of  pressure,  it  is  sufficient  to  define  the  amount  of  ozone  in  an 
atmospheric  path  with  a  single  parameter  proportional  to  the  total  number  of  mole¬ 
cules  in  the  path.  Figure  23  gives  the  ozone  amount  as  a  function  of  altitude  in 
(cm-atm)g^,p/km  for  each  of  the  five  atmospheric  models,  and  should  be  used  to 
determine  the  amount  of  ozone  In  a  horizontal  path.  Figure  24  provides  total  ozone 
amounts  between  the  top  of  the  atmosphere  and  the  indicated  altitudes  for  each 
model  atmosphere.  Ozone  amounts  for  vertical  layers  within  the  atmosphere  can 
be  found  by  subtraction  of  values  read  off  Figure  24  for  two  different  altitudes. 

Slant  path  amounts  are  found  by  multiplication  of  the  vertical  amount  by  sec  6  or 
M(0)  (see  Section  7).  For  a  discussion  of  the  unit,  (cm-atm)g^,p,  used  in  this 
section,  see  the  Appendix. 

The  amount  of  ozone  in  the  required  atmospheric  path  is  determined  from 
Figuren  23  and  24.  Then  the  atmospheric  transmittance  through  the  path  is  found 
as  h  function  of  wavelength  from  Figure  25  obtained  from  Inn  and  Tanaka  (1953). 

The  application  of  this  figure  is  identical  to  that  described  earlier  for  the  prediction 
charts  given  in  Section  5.4. 


6.  SCATTERED  SOLAR  RADIATION  AND  INFRARED  EMISSION 


6,1  Transmitted  and  Reflected  Sky  Radiance  in  Clear  and  Cloudy  Sky 

While  the  energy  taken  out  of  a  beam  of  radiation  by  absorption  contributes  to 
heating  of  the  air,  the  energy  scattered  by  molecules,  aerosols,  or  cloud  droplets 
will  be  redistributed  in  the  atmosphere.  We  will  only  consider  the  problem  of  sun¬ 
light  Incident  on  a  plane  parallel  atmosphere.  As  long  as  only  single  scattering  is 
Important  and  a  homogeneous  atmosphere  is  assumed  (aerosol/air  mixing  ratio 
constant  with  height),  the  angular  distribution  of  skylight  at  a  specific  wavelength 
can  be  obtained  from 
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(scattering  and  absorption  by  molecules  and  aerosols-see  Section  5.3 
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by  taking  the  logarithm  of  the  transmittance).  With  regard  to  molecular  . 
absorption  (km>,  Eq.  24  is  only  valid  for  monochromatic  radiation 
(see  Table  5  for  km  values) 


aa+ma  =t  +  t 
v  a  v  m  s,a  s,m 


F(4>)  3  angular  scattering  function  (molecular  +  aerosol)  per  unit  air  mass 
fg  a  +  ^m(4>)  ts  m  (see  Figure  26  for  ^a($)  and 

The  scattering  angle  ^  is  defined  by  cos  ^  =  cos  0  cos  01  -  sin  0  sin  ©'  cos  (A¥) 
whereA'I'  is  the  angular  azimuth  difference  between  the  sun  direction  and  the  line 
of  sight. 

However,  if  tM  or  tM1  becomes  larger  than  0.10,  higher  order  scattering  can 
no  longer  be  neglected;  this  means  that  illumination  of  the  scattering  volume  from 
the  sky  and  reflecting  earth's  surface  becomes  increasingly  important.  Higher 
order  scattering  becomes  dominant  for  tM  or  tM1  >0.5,  particularly  for  high* 
ground  albedo.  This  is  particularly  true  for  the  radiation  leaving  the  top 
of  the  atmosphere  (radiation  to  space).  The  general  radiative  transfer 
equation  has  been  solved  analytically  only  for  molecular  scattering, 
and  detailed  tables  of  sky  radiance,  polarization  (including  the  Stokes  parameters) 
have  been  prepared  by  Coulson  et  al  (1960).  Skyllghi  distributions  for  hazy  atmos¬ 
pheres  have  been  computed  by  de  Bary  et  al  (1965)  using  these  tables  and  Eq.  (24). 
However,  computers  have  made  possible  a  statistical  treatment  of  the  scattering 
problem  in  any  realistic  hazy  or  cloudy  atmosphere.  By  this  Monte  Carlo  method, 
several  investigators  have  computed  the  distribution  of  the  radiance  of  the  hazy 
atmosphere,  both  as  seen  from  the  ground  (Figure  27a)  and  from  space  (Figure  27b). 
The  assumed  aerosol  size  distribution  and  vertical  distribution  is  essentially  the 
same  as  in  Table  3.  Dotted  curves  of  the  radiance  (per  unit  solid  angle  and  unit 
Incident  solar  flux)  are  for  the  sun  in  the  zenith,  and  solid  lines  are  for  the  sun  at 
e  =  86.3°. 

In  Figure  27a,  the  downward  or  "transmitted"  radiance  generally  increases 
with  decreasing  wavelength  when  the  Bun  is  in  the  zenith.  Forward  scattering 
causes  the  radiance  to  peak  near  the  sun.  Near  the  horizon,  the  radiance  increases 
again,  except  at  short  wavelengths,  and  the  albedo  influence  is  large.  At  low  solar 
elevations,  downward  radiance  is  generally  much  smaller,  especially  near  the 
zenith  at  long  wavelengths,  and  the  albedo  influence  is  significant  only  at  short  wave' 
lengths.  Since  the  sky  radiance  values  have  been  averaged  over  intervals  of  0.  1 
in  cos  0  (that  is,  over  angular  Intervals  ranging  from  5°  to  25°),  the  radiance  near 
the  sun  actually  should  be  much  higher  than  shown. 

In  Figure  27(b),  upward  or  reflected  radiance  as  seen  from  outBide  of  the 
atmosphere  is  strongly  dependent  on  the  albedo  of  the  ground.  The  contribution  of 
reflected  ground  radiance  becomes  dominant,  especially  in  the  near  IR  and  in  the 
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nadir  direction  where  the  atmospheric  backscattered  flux  becomes  small.  This 
effect  will  be  less  pronounced  in  more  hazy  atmospheres. 

The  distribution  of  sky  radiance  in  the  ultraviolet  is  dominated  by  multiple 
scattering  and  below  0.  35  Mm  by  ozone  absorption.  At  high  solar  elevations,  down¬ 
ward  as  well  as  upward  radiances  at  0.  30  Mm  are  only  about  1  percent  of  their  res¬ 
pective  values  at  0.  40  Mm,  but  independent  of  albedo. 

In  the  case  of  the  radiance  from  clouds,  we  present  Figures  28a  and  28b  for 
a  rather  dense  nimbostratus  (Kattawar  and  Plass,  1968).  Aerosol  and  molecular 
scattering  have  been  neglected  in  this  model  because  they  would  have  little  effect. 
The  strong  downward  radiance  (Figure  28a)  for  high  sun  and  low  albedo  is  largest 
near  the  sun,  but  with  high  albedo  is  nearly  independent  of  the  angle  of  observation. 
At  low  sun,  the  much  smaller  cloud  radiance  is  largest  in  the  zenith.  The  radiance 
of  the  cloud  top  (Figure  28b)  again  is  very  large  at  high  sun.  At  low  sun  there  is 
litt  a  discrimination  between  a  cloudy  and  a  clear  sky  (see  Figure  27,  U0.7  jum). 

An  indication  of  the  variability  of  the  reflectance  or  spectral  albedo  of  some 
natural  objects  may  be  obtained  from  Figure  29.  The  data  presented  are  typical 
with  respect  to  the  wavelength  dependence,  but  absolute  reflectance  may  vary  con¬ 
siderably  for  individual  objects,  and  generally  is  higher  for  radiation  near  grazing 
incidence. 


6.2  Infrared  Emission  of  the  Atmosphere 

In  the  infrared  from  5  to  25  pm  and  beyond,  the  atmospheric  emission  in  a 
collimated  beam  can  be  determined  by  numerical  integration  of  the  equation  of 
radiative  transfer  to  give  the  upward  radiance,  IT  (v)\ 


IT  (v) 


B  (T,  T)  dT  +  B  (v,  T  )  t 

8  8 


(25) 


where  the  integral  represents  the  atmospheric  contribution  and  the  second  term  is 
the  contribution  from  the  ground  (or  lower  boundary)  modified  by  the  total  atmos¬ 
pheric  transmittance  between  the  boundary  and  the  observer;  r  is  the  transmit¬ 
tance  which  ranges  between  1. 0  at  the  level  of  observation  and  at  the  lower 
boundary  (for  example,  the  earth's  surface  or  cloud  top);  B  (v,  T)  is  the  Planck 
(bmck  body)  function  corresponding  to  frequency  v  and  the  temperature  T  of  an 
atmospheric  layer,  and  Tg  is  the  ground  temperature. 

To  obtain  the  downward  emission  from  a  clear  sky,  the  second  term  in  Eq,  (25) 
should  be  omitted  and  is  then  the  transmittance  of  the  entire  atmosphere.  In  a 
cloudy  sky  the  subscript,  g,  refers  to  a  cloud  base. 

The  effects  of  aerosols  and  of  molecular  scattering  have  been  neglected  and 
the  molecular  absorption  charts  described  in  Section  ft.  2  have  been  applied  to  the 
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water  vapor,  ozone,  and  uniformly  mixed  gases  corresponding  to  the  Midlatitude 
Winter  model.  Equation  (25)  was  then  applied  using  the  derived  transmittance 
values  and  performing  the  numerical  integration  according  to  a  numerical  analogue 
of  Eq.  (25). 

The  results  are  presented  in  Figure  30  for  both  downward  and  upward  beam 
radiation.  The  downward  radiation  case  assumes  the  observer  to  be  at  the  ground 
looking  up  at  the  zenith.  The  upward  case  assumes  an  observation  point  at  the  top 
of  the  atmosphere  viewing  the  nadir  (that  is,  looking  down). 

It  should  be  noted  that  in  the  upward  case  the  major  absorption  features  (6.  3 /4m 
H20,  9.  6  pm  O3,  and  15  pm  C02)  appear  as  relative  minima,  whereas  in  the  down¬ 
ward  case  these  same  regions  appear  as  relative  maxima.  The  surface  tempera¬ 
ture,  Tg,  was  taken  to  be  the  same  as  the  temperature  of  the  lowest  atmospheric 
level  in  the  model  (272K).  The  upper  curve  in  Figure  30  corresponds  to  a  black- 
body  at  this  boundary  temperature. 

For  wavelengths  shorter  than  5  pm,  the  problem  of  upward  and  downward 
radiance  is  complicated  by  molecular  and  aerosol  scattering  and  is,  therefore,  a 
function  of  the  solar  zenith  angle.  Computations  of  radiation  between  3  and  5 
micrometers  is  particularly  complicated.  Some  examples  of  measured  radiation 
in  the  1  to  5  pm  region  are  included  in  Figure  31. 


7.  REFRACTIVE  EFFECTS  IN  THE  ATMOSPHERE 

Previously  (Sections  3  and  5)  it  was  indicated  that  slant  path  length  could  be 
obtained  from  vertical  path  length  through  multiplication  by  sec  9  for  6  <  80°.  For 
9  &80°,  the  vertical  path  length  (or  other  concentration  parameter),  AL,  may  be 
converted  to  slant  path  length  by 

AL  (slant  path)  (Vertical)  X  M(0)  (26) 

where  M(0)  is  the  relative  optical  path,  defined  as  the  ratio  of  slant  path  length  to 
vertical  path  length  between  two  levels.  The  term  M(9)  always  equals  sec  9  for  a 
short  path  near  the  observer,  and  for  any  path  if  6  <  80°.  For  a  long  slant  path, 
k.wY.over,  M(0)  becomes  smaller  than  sec  6  because  of  the  curvature  of  the  atmos¬ 
phere,  which  is  partly  compensated  by  refraction  of  the  beam.  This  implies  that 
the  vertical  distribution  of  a  scattering  or  absorbing  atmospheric  constituent  also 
Influences  the  optical  path,  and  the  slant  path  length  may  have  to  be  derived  from 
Eq.  (26)  if  high  accuracy  is  warranted. 

If  curvature  and  refraction  are  considered,  the  general  equation  for  the  optical 
path  between  altitudes  iQ  and  z  ^  of  a  constituent  of  density  p  at  z  and  scale  height 


# 
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H  =  (Zj  -  Zq)  /loge  £  p(zQ)/p(Zj)J  may  be  written 


1-1  fz 

f  /r+znf 

/  n„  \?1 

i 

b  -hfl 

(nM  sin  ®)  J 

-1/2 


p(z)  dz  (27) 


where  n  is  the  refractive  index  of  air  and  the  (o)  subscript  refers  to  the  lower 
reference  level,  and  r  is  the  radius  of  the  earth. 

For  the  molecular  atmosphere  and  uniformly  mixed  atmospheric  constituents, 
the  air  mass  relative  to  the  vertical  air  mass  of  the  entire  atmosphere  measured 
from  the  ground,  M(0),  as  computed  by  Fasten  (1967)  is  presented  in  Figure  32. 
These  data  are  based  on  the  ARDC  Model  Atmosphere  of  Minzner  et  al  (1959)  and 
a  refractive  index  of  air  corresponding  to  a  wavelength  of  0.  7  pm.  Earlier  values 
of  M(0)  of  which  the  data  by  Bemporad  (1907)  are  most  often  used,  are  in  fair 
agreement.  Fasten  (1967)  has  shown  that  the  average  influence  of  latitude  and 
season  onJVl(0)  is  less  than  2  percent  but  abnormal  tropospheric  density  profiles 
may  have  larger  effects.  Figure  32  compares  M(0)  values  for  dry  air  with  M(0) 
values  for  an  atmosphere  containing  water  vapor  with  a  distribution  according  to 
Schnaidt  (1938).  Values  of  sec  0  are  included  for  comparison 

The  angular  deviation  between  the  actual  and  apparent  solar  zenith  angles,  A 0. 
due  to  refractive  effects  is  presented  in  Table  6. 


Table  6.  Difference  between  Actual  and  Apparent  Zenith  Angle  as  a  Function 
of  Apparent  Zenith  Angle 


%p(DEG) 

75 
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83 
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0.  124 

0.  167 

0.  199 

0.248 

0app  (DEG) 

88 

88.3 

89 

89.5 

00 

60  (DEG) 

0.311 

0.358 

0.417 

0.497 

0.  570 
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For  thin  atmospheric  layers  (10  km  or  less)  at  high  altitudes  (where  refractive 
effects  can  be  neglected).  Table  7  provides  the  geometric  portion  (that  is,  the  ef¬ 
fect  is  due  to  curvature  only)  of  the  slant  path  correction.  For  example,  the  geo¬ 
metric  path  through  a  thin  layer  at  30  km  viewed  at  a  zenith  angle  of  00°  is  10.  36 
times  the  geometric  path  in  the  vertical.  However,  the  optical  path  through  a 
uniformly  mixed  atmospheric  constituent  at  the  same  viewing  angle  would  be  36.  3 
(see  Figure  32).  This  situation  may  occur  in  the  treatment  of  transmittance  through 
the  ozone  layer  or  stratospheric  dust  layers. 


Table  7.  Relative  Optical  Path  of  a  Thin,  High  Altitude  Layer  as  a 
Function  of  Zenith  Angle  9 


Height  of  Layer  (km)  Above  Surface 
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10.  30 

8.  94 

8.  01 

7.  33 

90 
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12.  66 
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8.  97 

8.03 

7.  34 

8.  EFFECT  OF  Cl, 01  OS  AMI)  FOCI  ON  ATMOSPHERIC  RADIATION 

As  in  the  case  of  aerosol  extinction  (Section  3.  2)  when  a  beam  of  radiation  is 
intercepted  by  a  cloud  or  fog  layer,  some  of  the  incident  radiation  may  be  absorbed 
and  some  may  be  scattered  out  of  this  beam  with  both  effects  thereby  causing  a  net 
attenuation  of  the  original  radiation.  The  magnitude  of  the  attenuation  depends  on 
the  number,  size  distribution  and  refractive  index  of  the  suspended  droplets  and 
on  the  thickness  of  the  layer.  Since  the  radiation  may  be  scattered  in  all  direc¬ 
tions,  including  those  close  to  that  of  the  original  beam,  the  resulting  attenuation 
may  also  depend  on  the  acceptance  angle  of  the  receiver. 

Natural  fogs  and  low  level  clouds  are  composed  of  spherical  water  droplets, 
the  refractive  properties  of  which  have  been  fairly  well  documented  in  the  spectral 
region  of  interest  here  (Kislovskii,  1959;  and  irvine  and  Pollack,  1966).  However, 
in  the  case  of  cirrus  clouds  the  particles  consist  of  ice  crystals  of  various  shapes 
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(hexagonal  or  cylindrical  crystals)  which,  when  mutually  oriented,  can  give  rise 
to  anomalous  scattering  and  polarization  phenomena  when  radiation  is  propagated 
through  them. 

In  this  section  we  will  restrict  outselves  to  the  propagation  of  a  collimated 
beam  of  radiation.  Generally,  fogs  and  clouds  show  only  a  small  increase  in  trans¬ 
mittance  with  wavelengths  from  the  visible  to  about  6  Mm.  Beyond  6  urn,  a  greater 
increase  is  observed.  However,  certain  infrequently  occurring  types  of  fog 
characterized  by  a  sharply  peaked  distribution  of  small  droplets  (this  phenomenon 
is  also  found  near  the  base  of  stratus  and  cumulus  clouds),  show  a  significant  in¬ 
crease  in  infrared  transmittance. 

These  general  effects  can  be  seen  from  Figure  33,  which  shows  the  wavelength 
variation  of  the  extinction  coefficient  for  a  drop  size  distribution  representative  of 
fair  weather  cumulus  cloud  after  Deirmendjian  (1964  and  19C9).  The  model  drop 
size  distribution  (shown  in  Figure  34),  which  is  given  by 

n(r)  =  2.  373  r6  e  **  cm  ^  Mm  *,  (28) 

has  a  maximum  corresponding  to  a  radius  of  4  Mm  and  i3  normalized  to  100  drops 
cm  "3. 

Figure  33  was  calculated  and  based  on  single  scattering  theory,  using  the 
complex  refractive  index  values  for  water  given  by  Kislovskii  (1959)  and  Eq.  (28). 

It  will  be  seen  from  Figure  33  that  the  attenuation  due  to  this  fairly  thin  model 
cloud  is  rather  severe;  for  example,  for  a  100m  path  length  the  transmittance  is 
approximately  0.  3  and  0.  15  at  10  Mm  and  1  r  >  respectively. 

Some  mention  should  be  made,  at  this  point,  of  the  validity  of  using  single 
scattering  theory  to  account  for  the  attenuation  of  radiation  in  cloud. 

It  has  generally  been  assumed  that  calculations  of  attenuation  based  on  single 
scattering  theory  are  valid  for  an  optical  thickness  of  less  than  0.  03  (van  de  Hulsts 
criterion)  and  that  for  values  of  optical  thickness  in  excess  of  unity,  secondary  and 
multiple  scattering  effects  become  important.  However,  this  applies  to  the  propa¬ 
gation  of  diffuse  radiation.  For  a  narrow  beam  of  collimated  radiation  (or  laser 
b  .am  propagation)  it  has  been  found  (Zuev  et  al,  1967)  that  single  scattering  theory 
can  account  for  the  observed  attenuation  for  values  of  opiical  thickness  of  up  to  2  5. 

In  this  context  optical  thickness  is  defined  as  the  product  of  the  extinction  coeffi¬ 
cient  and  path  length. 
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Figure  3.  Equivalent  Sea  Level  Path  Length  of  Water  Vapor  as  a 
Function  of  Altitude  for  Vertical  Atmospheric  Paths 
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Figure  8.  Equivalent  Sea  Level  Path  Length  of  Nitrogen  as  a  Function 
of  Altitude  for  Horizontal  Atmospheric  Paths 
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Ft<r4re  9.  Equivalent  Sea  Level  Path  Length  of  Nitrogen 
of  Altitude  for  Vertical  Atmospheric  Paths 
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Figure  12.  Equivalent  Sea  Level  Path  Length  for  Molecular  Scattering 
aa  a  Function  of  Altitude  for  Horizontal  Atmoapherlc  Pa  the 
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Figure  13.  Equivalent  Sea  Level  Path  Length  for  Molecular  Scattering 
as  a  Function  of  Altitude  for  Vertical  Atmospheric  Paths 
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Figure  14.  Equivalent  Sea  Level  Path  Length  for  Aeroeol  Extinction 
aa  a  Function  of  Altitude  for  Horliontal  Atmospheric  Paths 
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Figure  16b.  Prediction  Chart  for  Water  Vapor  Transmittance  (4-26  frni) 
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Figure  18b.  Prediction  Chart  for  Ozone  Transmittance  (8-18  (im) 
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Prediction  Chart  for  Transmittance  Due  to  Molecular  Scattering 
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Figure  23.  Equivalent  Sea  Level  Path  Length  for  Ozone  (0.25-0.  75  fan)  as  a  Function 
for  Horizontal  Atmospheric  Paths 
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Figure  29.  Typical  Reflectance  of  Water  Surface,  Snow,  Dry  Soil  and  Vegetation  (see 
Krinov,  1953;  Penndorf,  1956;  and  Earing  and  Smith,  1965) 
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Figure  31.  Spectral  Radiance  of  a  Clear  Zenith  Sky  Showing  the  Effect  of  the 
Position  of  the  Sun  on  the  Radiation  Scattered  Below  3  /um  and  the  Effect  of 
Thermal  Emission  at  Longer  Wavelengths.  Measured  at  Cocoa  Beach, 
Florida  (see  Bell  et  al,  1960) 
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Figure  32.  Comparison  of  the  Relative  Air  Mass  [M(0)]Computed  by  Fasten  (1967) 
as  a  Function  of  Apparent  Zenith  Angle,  9,  Based  on  the  A RllCLModel  Atmosphere 
(Mlnzner  et  al,  1959)  With  a  Similar  Quantity  [Mw(0)J  Based  on  a  vertical  Water 
Vapor  Distribution  According  to  Schnaldt  (1938)  and  for  sec  9  (Note  that  for  0<6O°, 
m(9)  -  M w(0)  =  sec  9.) 
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Attenuation  as  a  Function  of  Wavelength  Due  to  Fair-weather  Cumulus  Cloud 
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Appendix 


1.  INTRODUCTION 

Throughout  thin  report  we  have  used  the  term  "frequency"  rather  loosely  to 
signify  spectral  variation.  In  general,  frequency  is  assumed  to  he  synonymous 
with  wavenumber,  v  (cm-1),  which  is  given  by  v  ■  104A  where  the  wavelength  A 
is  ^iven  in  units  of  /um  (that  is,  10'®  meter). 

The  subject  of  units  used  in  atmospheric  physics  is  rather  confusing  to  the 
beginner.  This  is  iwrtlcularly  true  where  units  of  absorber  concentration  or 
amount  are  concerned.  In  this  Appendix  we  will  show  the  reader  how  to  derive  the 
absorber  concentrations  for  use  with  the  prediction  charts  given  in  Section  5,  on 
the  assumption  that  one  has  the  relevant  meteorological  data  available.  We  will 
then  discuss  the  units  of  absorber  concentration  and  indicate  how  to  make  the  ap¬ 
propriate  conversions. 


2.  ABSORBER  CONCENTRATION 

For  most  gases  absorber  concentration  denoted  here  by  AL  is  generally  mea¬ 
sured  by  the  equivalent  length  of  pure  absorber  in  the  path  under  the  existing  condi¬ 
tions  of  pressure,  P,  and  temperature,  T.  The  unit,  centimeter  atmospheres, 

(cm -atm)  is  given  by  AL  ■  cRP  where  c  is  the  fractional  concentration  of  the  ab¬ 
sorber  (see  Table  2),  K  is  the  path  length  (om)  and  the  P  is  the  pressure  in  atmos¬ 
pheres.  More  often  AL0  Is  quoted  as  the  absorber  concentration  reduced  to  STP 


7  ^  ' 


V 


A2 


conditions,  which  is  a  more  meaningful  quantity  where  vertical  or  slant  atmospheric 
paths  are  concerned.  Thus,  for  a  horizontal  path  at  altitude,  z,  we  have 


ALQ  ■  cR 


(cm -a  tin  Ig'pp 


(Al) 


where  TQ  ■  273K,  and  PQ  ■  1  atm.  If  c  is  given  in  part?  per  million  (ppm),  (see 
Table  2),  and  R  is  given  in  km,  the  following  expression  can  be  used. 


*Lo 


(cm-atm)g^,p  =10  *cR 


(A2) 


For  a  uniformly  mixed  gas  of  fractional  concentration  by  volume,  c,  the  number  of 
centimeter  atmospheres  at  STP  in  a  slant  path  from  any  altitude  z  to  the  limit  of 
the  atmosphere  is  given  by 


AL  (6,  z  — •  oo  ) 


STP 


c  H  (z)  sec  9 


P(z) 

P„ 


(to)- 


cH(z)  sec  9  p( z)  /  p 


(A3) 


where  H(z)  is  called  the  scale  height  of  the  atmosphere  above  altitude  z,  9  is  the 
zenith  angle  and  p(z)  and  pQ  are  the  atmospheric  density  at  altitude  z  and  STP, 
respectively.  For  9  >  70°,  see  Section  7.  From  the  hydrostatic  equation,  H(z) 
can  be  defined  as  follows: 

P(z)  =  gp  (z)  H(z)  (A4) 

where  P(z)  is  the  pressure  at  altitude  z,  and  g  is  the  acceleration  due  to  gravity 
(981  cm  sec'2).  Thus,  for  a  vertical  column  of  air  at  STP,  H(o)  ■  7.  99  x  105  cm 
(that  is,  approximately  8  km).  If  the  scale  height  in  cm  is  required  at  another 
altitude,  z,  it  can  be  determined  from 

H(z)  ■  1.  02  x  106  x  [P(z)/p(z)J  (A5) 

where  P(z)  is  given  in  millibars  and  p(z)  is  given  in  gm/rn2  (see  Table  1).  Note 
from  Eq.  (A5)  that  the  quantity  H(z)  p(z)  is  directly  proportional  to  P(z)  and  is 
independent  of  temperature.  inserting  H(z)  p(z)  into  Eq.  (A3),  we  have 

ALq  (9,  z  -  « )STp  •  c  x  1.  02  x  10®  P(z)  sec  9  /pQ  (A6) 
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and  for  the  amount  of  absorber  between  levels  Zj  and  z 2  (see  Table  2): 

^LQ  (6,  Zj  —  z2^STP  “  c  *  L  02  x  1q6  I  p(zi>  ‘  p<z2^  sec  9/P o'  (A7) 

For  condensible  vapors  (for  example,  water  vapor)  it  is  more  convenient  to 
measure  the  total  length  of  liquid  (in  ems)  which  may  be  precipitated  out  of  the  path 
(per  unit  area).  The  number  of  precipitable  centimeters  (abbreviated  pr  cm)  of 
water  vapor  in  a  given  path  of  length  R(cm)  is  given  by,  w  ■  p(T)  Wr  where  p(T)  is 
the  saturation  vapor  density  of  water  at  temperature  T  (gm  cm*2)  and^  is  the 
fractional  relative  humidity.  It  will  be  noted  that  the  units  of  pr  cm  and  gm  cm'2 
are  equivalent. 

Figure  A1  will  enable  the  reader  to  calculate  the  number  of  pr  ems  of  water 
vapor  per  km  path  length,  for  a  given  temperature  and  relative  humidity.  For 
temperatures  below  -40°C  the  reader  should  use  density  values  of  water  vapor 
measured  oyer  ice  for  determining  water  vapor  concentrations  for  atmospheric 
paths  (see  Smithsonian  Meteorological  Tables.  1963). 
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In  Table  1,  water  vapor  densities  are  given  in  units  gm  m  °  as  a  function  of 
altitude  for  the  5  model  atmospheres.  It  may  be  useful  to  remember  that  for  water 
vapor, 

I  gm  m'3  ■  0.  1  gm  cm'2/km  (or  0.  1  pr  cm/km). 


3.  UNIT  CONVERSION  FACTORS 

According  to  Avogadros  hypothesis  the  molecular  weight,  M,  in  gms. ,  or  any 
gas  occupies  22.4  liters  at  STF.  Hence  one  cubic  centimeter  of  gas  at  STP  weighs 
(M/2.  24,  104)  gms.  Since  1  (cm -atm )g  1?,  is  equivalent  to  a  length  of  1  cm  of  gas 
at  STP  per  cm2,  we  can  write  the  following: 

1  cm -atm  of  gas  ■  - —  — —  gm  cm-2  (A8) 

2.24  x  104 

where  M  Is  the  molecular  weight  of  the  gas  (see  Table  2). 

Also  since  1  cm3  of  gas  at  STP  contains  1/2.  24x10*  moles  and  one  mole  of  gas 
contains  the  same  number  (&)  of  molecules  (where  (PL  is  Avogadros  Number)  we 
can  write  that  1  cm3  of  gas  at  STP  contains  6.  06x  1023/2. 24x  10*  •  2.  60x  101® 
molecules  or 


1  (cm -atm). 


2.89x  10 18  molecules/cm2. 


(A9) 


A4 


Note  that  Eq.  (A9)  la  valid  independent  of  the  gas.  Since  for  water  vapor  one  can 
write  1  gm  cm’2  in  terms  of  (cm-atm)gTp,  using  Eq.  (A8)  it  will  also  be  seen  that 
for  water  vapor: 

1  pr  cm  (HgO)  ■  3.  34x  1022  molecules/cm2.  (A10) 

Thus  the  unit,  molecule/cm2,  is  independent  of  the  nature  of  the  absorbing  gas  and 
basic  to  all  gases.  It  is  hoped  that  this  unit  will  become  more  widely  used  and 
thereby  help  to  unify  absorber  concentration  units  for  all  the  atmospheric  consti¬ 
tuents. 
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Figure  Al.  Water  Vapor  Concentration  per  Kilometer  Path  Length  as  a  Function  of  Temperature 
and  Relative  Humidity 


